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output  of  the  beamformer,  signal-to-noise  ratio,  and  probability  of  detection.  Reverberation  plots  show  surface, 
bottom,  and  volume  components  as  well  as  total  reverberation  and  competing  noise  level. 
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1000  kyd;  a  maximum  of  500  target  ranges  are  flowed  per  run.  Source  frequency  is  allowed  to  be  between  25 
Hz  and  25  kHz.  Pulse  length  may  be  be  tween  (f0^7  and  100  s. 

The  sound-speed  profile  is  represented  by  curvilinear  segments  to  eliminate  false  caustics.  The  sound-speed 
profile  is  assumed  to  be  constant  with  range  and  the  ocean  bottom  horizontal. 

The  propagation  loss  models  are-fljfray  theory  supplemented  with  caustic  corrections  derived  by 
Brekhovskikh,  and42)‘empirical  equations  based  on  the  AMOS  data  and  corrected  for  low  frequencies  using 
normal  mode  theory.  Beam  patterns  for  both  transmitter  and  receiver  are  incorporated  into  propagation  loss. 
The  Hall-Watson  model  is  used  for  absorption  loss. 

Surface  backscattering  strength  is  a  combination  of  the  Chapman-Harris  equations,  Eckart’s  equations,  and 
Richter’s  data.  Bottom  backscattering  strength  is  derived  from  Lambert’s  law  and  Schmidt’s  data.  Volume 
backscattering  strength  is  represented  as  the  column  strength.  ^ - - 

Reverberation  is  averaged  over  the  pulse  length.  Contributions  from  previous  pings,  if  any,  are  included  in  the 
reverberation  level  for  the  current  ping.  Reverberation  is  corrected  for  doppler  gain  due  to  target  motion. 
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for  various  assumptions  of  signal  distortion  and  detector  characteristics  are  available. 

The  computer  program  uses  10  -  300  s  of  cpu  time  per  run  on  the  UNIVAC  1 1 10  computer  and  requires 
57,000  words  of  core  storage. 
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SUMMARY 


OBJECTIVE 

Extend  the  capabilities  of  LORA  (Ref.  1),  a  computer  model  for  predicting  the  per¬ 
formance  of  active  sonars,  to  simulate  active-sonar  surveillance  scenarios.  Specifically,  include 
the  following  additional  capabilities: 

1.  Frequency  regime  shall  be  25  Hz  to  25  kHz. 

2.  Broad  bandwidths  and  long  pulse  lengths  for  FM  and  PRN  signals  shall  be  allowed. 

3.  Target  range  shall  be  increased  to  a  1000-kyd  maximum. 

4.  Source  and  receiver  shall  not  be  restricted  to  the  same  depth. 

5.  Output  plots  shall  be  made  available. 

COMPUTER  PROGRAM  DESCRIPTION 

The  new  computer  program,  LIRA,  contains  many  of  the  capabilities  of  its  pre¬ 
decessor,  LORA.  The  capabilities  of  LIRA,  after  incorporating  the  new  additions  to  the 
model,  include: 

1 .  Active-sonar  performance  predictions  can  be  obtained  for  target  ranges  extend¬ 
ing  from  .001  kyd  to  1000  kyd. 

2.  Source  frequencies  can  be  from  25  Hz  to  25  kHz. 

3.  The  source  and  receiver  are  not  restricted  to  the  same  depth.  (The  target  can 
be  at  still  another  depth). 

4.  The  scenario  of  a  receiving  array  mounted  on  a  steeply  sloping  shelf  is  modeled. 

5.  Curvilinear  techniques  are  used  with  the  sound-speed  profile  to  eliminate  false 
caustics. 

6.  Propagation  losses  calculated  using  ray  theory  are  supplemented  with  caustic 
corrections  derived  from  wave  theory. 

7.  The  spatial  positions  of  caustics  are  determined  by  iterative  methods. 

8.  Beam  patterns  for  both  transmitter  and  receiver  are  allowed. 

9.  The  surface  duct  propagation  loss  model  is  based  on  the  AMOS  equations 
modified  for  low  frequencies  using  normal  mode  theory. 

10.  The  Hall-Watson  model  is  used  for  absorption  loss. 

1 1 .  The  surface  backscattering  strength  is  derived  from  the  Chapman-Harris 
equation,  Eckart’s  equations,  and  Richter’s  data  (the  user  may  also  enter  his  own  model). 

12.  Bottom  backscattering  strength  is  derived  from  Lambert’s  law  and  Schmidt’s 
data  (the  user  may  also  enter  his  own  model). 

13.  Volume  backscattering  strength  is  represented  as  the  water  column  back- 
scattering  strength. 
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14.  Reverberation  is  averaged  over  the  pulse  length. 

1 5.  A  maximum  pulse  length  of  100  s  is  allowed. 

16.  Reverberation  contributions  from  previous  pings,  if  any,  are  included  in  the 
reverberation  level  for  the  current  ping. 

17.  Doppler  gain  against  reverberation  can  be  calculated  by  the  program  or 
entered  as  an  input. 

18.  CW,  FM,  and  PRN  signals  are  allowed. 

19.  Reverberation  and  noise  are  corrected  for  processing  gain. 

20.  Either  recognition  differential  (at  the  output  of  the  beamformer)  or  detection 
threshold  (at  the  output  of  the  signal  processor)  may  be  used  to  represent  the  detection- 
decision  process. 

21 .  Signal  excess  can  be  calculated  using  the  propagation  loss  to  the  target  for 
either  the  incoherent  sum  of  ray  arrivals  or  the  largest  arrival. 

22.  Five  probability-of-detection  models  are  available  for  various  assumptions  of 
signal  distortion  and  detector  characteristics. 

23.  Plots  are  available  for  propagation  loss,  reverberation,  and  signal  excess.  Rever¬ 
beration  plots  display  surface,  bottom,  and  volume  components  as  well  as  total  reverbera¬ 
tion  and  competing  noise  level.  Propagation  loss  and  signal  excess  plots  each  contain  two 
curves,  one  for  the  incoherent  sum  of  ray  arrivals  and  one  for  the  largest  arrival. 

LIRA  has  the  following  limitations: 

1.  The  bottom  is  assumed  to  be  horizontal. 

2.  The  sound-speed  profile  cannot  vary  with  horizontal  range. 

3.  The  source  and  receiver  cannot  be  displaced  in  either  range  or  azimuth. 

The  LIRA  computer  program  has  been  executed  thousands  of  times  and  is  able  to 
perform  many  successive  runs  without  failure.  Computer  costs  are  low,  from  $2  to  $1  5 
per  run  (10  to  300  s  of  cpu  time  for  the  UNI  VAC  1110  computer).  The  average  cost 
for  surveillance  predictions  is  S5  with  plots  extra  (about  S4  per  plot).  Computer  core 
storage  required  is  57,000  words.  The  program  is  segmented  to  allow  for  plotting. 

RECOMMENDATIONS 

Caustics  in  deep  submerged  ducts  give  rise  to  excessively  low  propagation 
losses  for  certain  geometries.  Either  a  model  to  describe  leakage  of  sound  from  the 
duct  should  be  developed  and  validated  or  another  model  currently  available  should 
be  chosen  for  the  incorporation  into  LIRA. 

Separation  of  source  and  receiver  in  range  should  be  modeled  as  the  next  step 
in  the  development  of  a  bistatic  model. 

Plots  of  normalized  reverberation  vs  time  should  be  provided  so  that  comparisons 
with  other  reverberation  models  can  be  made  directly. 

The  fields  at  horizontal  caustics  should  be  investigated. 
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INTRODUCTION 


LIRA,  Low-frequency  Intermediate  Range  Active-sonar  performance  prediction 
program  is  a  revision  and  an  extension  of  LORA,  LOng-Range  Active-sonar  performance 
prediction  program  (Ref.  1 ).  The  LORA  computer  program  was  created  to  predict  the 
performance  of  active-sonar  systems  for  frequencies  above  1  kHz  and  for  ranges  not 
greater  than  five  convergence  zones.  For  active-sonar  systems  using  this  frequency 
regime,  three  convergence  zones  is  long  range. 

Stimulated  by  active-source  technologists,  the  Navy  in  1977  renewed  its  interest 
in  mid-  to  low-frequency  active-sonar  surveillance.  Two  modeling  requirements  arose: 

1 .  Determine  the  optimum  design  parameters  of  a  new-technology  active  source, 
i.e.,  specify  source  level,  bandwidth,  and  frequency.  Highest  priority  must  be  given  to 
using  current  receiver  assets. 

2.  Make  performance  predictions  for  combinations  of  operating  environments 
and  system  parameters. 

LIRA  was  developed  to  play  a  role  in  satisfying  both  of  these  requirements. 

LORA  was  inadequate  for  the  following  reasons: 

1 .  LORA’s  active-sonar  performance  prediction  capability  is  limited  to  five 
convergence  zones  and  five  bottom  bounces.  The  maximum  range  for  accurate  predic¬ 
tions  is  therefore  limited  to  75  kyd  in  upward-refracting  areas,  such  as  the  Mediterranean, 
and  areas  characterized  by  a  submerged  double  duct. 

2.  LORA  does  not  have  caustic  corrections.  Below  i  kHz  the  acoustic  fields  from 
caustics  extend  noticeably  into  the  shadow  zones  predicted  by  classical  ray  theory. 

Below  100  Hz  these  fields  dominate  a  curve  of  propagation  loss  vs  range. 

3.  LORA  does  not  have  the  capability  to  generate  output  plots. 

4.  LORA’s  signal  processing  model  is  designed  for  a  short  CW  pulse.  The 
possibility  of  using  a  PRN  waveform  with  a  long  pulse  for  active  surveillance  requires 
averaging  the  reverberation  over  pulse  length  and  correcting  both  reverberation  and  noise 
for  the  time-correlation  gain. 

• 

Because  of  these  deficiencies  with  respect  to  the  modeling  requirements,  LORA 
did  not  impact  initial  system-design  studies.  However,  the  decision  was  made  to  extend 
LORA  to  low  frequencies,  longer  ranges,  and  longer  pulse  lengths.  A  capability  to  model 
bistatic  sonar  (separation  of  source  and  receiver  in  three  dimensions)  was  also  requested, 
but  constraints  of  time  and  money  allowed  only  the  modeling  of  depth  separation  of 
source  and  receiver. 

The  body  of  this  report  describes  the  LIRA  model  in  those  aspects  in  which  it 
differs  from  its  predecessor,  LORA,  which  is  described  in  Ref.  1.  Some  of  the  capabil¬ 
ities  of  the  LIRA  program  that  are  mentioned  in  the  Summary  are  not  described  in  the 
text  of  this  report  because  they  are  described  in  Ref.  1 .  The  topics  that  are  included 
here  are  restructured  ray  tracing,  caustics,  pulse-averaged  reverberation,  and  source- 
receiver  depth  separation. 
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The  first  section  describes  the  revised  propagation  loss  model  used  in  LIRA.  Basic 
definitions  are  given  for  clarity.  The  logical  order  of  steps  used  in  the  program  to  calculate 
propagation  loss  by  merging  ray  theory  with  caustic  corrections  derived  from  wave  theory 
is  presented.  The  properties  of  caustics  are  described  and  equations  are  given  for  the  pro¬ 
pagation  loss  at  and  near  caustics.  LIRA’S  methods  to  isolate  caustics  and  to  determine 
iteratively  the  caustic  parameters  are  described.  Tne  next  section  describes  the  reverbera¬ 
tion  model  as  it  differs  from  that  described  in  Ref.  1 .  The  addition  of  long  pulse  lengths 
to  LIRA’S  list  of  options  requires  that  the  reverberation  be  averaged  over  the  time  interval 
(pulse  length)  that  the  signal  is  being  received.  Pulse-averaged  reverberation  is  described 
for  CW.  FM,  and  PRN  signals.  The  final  section  deals  with  the  separation  of  source  and 
receiver  in  depth.  It  is  shown  that  propagation  loss  for  all  source-target-receiver  paths 
is  equivalent  to  the  average  of  the  two-way  propagation  losses  from  source  to  target  and 
from  receiver  to  target.  Reverberation  is  similarly  represented,  and  the  errors  introduced 
by  using  monostatic  backscattering  strengths  instead  of  bistatic  backscattering  strengths 
are  estimated.  A  summary  of  model  and  software  improvements  for  LIRA  is  given. 

Results  of  comparisons  with  other  models  (Appendix  D)  are  described.  Recommendations 
for  further  work  are  given. 

Appendix  A  presents  the  logical  structure  of  the  LIRA  program  to  aid  the  user  in 
streamlining  his  inputs.  Appendix  B  gives  descriptions  of  inputs  for  the  LIRA  program, 
i.e.,  definitions  of  parameters,  acceptable  limits  of  parameters,  typical  values,  etc. 

Appendix  C  describes  the  input  format  in  detail;  examples  are  given.  In  Appendix  D 
LIRA'S  propagation  loss  calculations  are  compared  with  those  from  PLRAY,  FACT,  and 
AP2.  Appendix  E  gives  an  example  of  using  LIRA  for  predicting  performance  in  a 
scenario  in  which  the  source  is  located  above  a  receiving  array  lying  on  a  sloping  shelf. 
Detection  ranges  are  printed  and  output  plots  are  given  for  propagation  loss,  reverberation, 
and  signal  excess. 


REVISED  PROPAGATION  LOSS  MODEL 


The  propagation  loss  model  used  in  LORA  has  been  revised  to  incorporate  caustic 
corrections.  This  section  explains  the  approach  used  in  the  LIRA  program  for  merging 
ray  theory  with  the  expressions  for  caustic  fields.  First,  some  basic  definitions  are  given 
to  clarify  terminology  used  throughout  the  report.  Then,  the  procedure  used  by  the  LIRA 
program  to  trace  rays,  isolate  caustics,  and  merge  ray  and  caustic  propagation  losses  is 
presented  in  a  step-by-step  format.  Following  this,  the  discussion  of  caustics  is  greatly  ampli¬ 
fied.  Equations  are  given  for  the  caustic  fields.  The  iterative  procedure  to  determine  caustic 
parameters  exactly  is  described  in  detail. 

BASIC  DEFINITIONS 

Throughout  this  report,  reference  is  made  to  ray  paths,  ray  segments,  ray  para¬ 
meters,  ray  types,  and  ray  families.  Definitions  of  these  terms  will  facilitate  the  dis¬ 
cussion  which  follows.  The  definitions  are  given  in  order  to  avoid  confusion  with  alterna¬ 
tive  usage  by  other  authors.  Because  ray  tracing  is  described  in  detail  in  the  companion 
document  describing  LORA  (Ref.  1 ,  pp.  25-34),  the  definitions  given  here  are  only 
supplementary. 

A  ray  path  is  the  path  along  which  sound  travels  through  the  water  from  one 
point  of  interest  to  another.  The  ray  path  may  contain  turning  points,  i.e.,  reversals  of 
direction  in  depth,  either  refractive  or  reflective.  The  depth  at  which  a  turning  point 
occurs  is  called  the  vertex  depth. 

A  ray  segment  is  a  part  of  a  ray  path  connecting  two  depths  of  interest.  Ray  paths 
connecting  two  depths  are  constructed  by  assembling  ray  segments.  A  ray  segment  has  no 
complete  turning  points,  although  an  end  of  a  segment  may  be  at  a  turning  point.  In  LIRA 
three  basic  ray  segments  are  used:  upper  turning  point  to  source  depth,  upper  turning 
point  to  target  (or  backscattering)  depth,  and  upper  turning  point  to  lower  turning  point. 

The  last  segment  is  one-half  of  a  ray  cycle. 

Ray  parameters  comprise  the  range  and  its  derivatives,  travel  time,  and  absorption 
loss.  Ray  paths  and  ray  segments  may  be  thought  of  as  vectors  with  ray  parameters  as 
components.  A  ray  path  is  the  vector  sum  of  ray  segments  appropriately  chosen  so  that 
the  ray  path  connects  two  depths  of  interest. 

A  ray  type  is  one  of  four  possible  ways  that  ray  paths  can  be  constructed  using 
combinations  of  the  same  ray  segments  to  connect  two  depths.  The  four  ray  paths  so 
constructed  have  their  end  points  at  the  same  depths,  but  the  ray  parameters  are,  in 
general,  different.  The  angles  of  the  ray  paths  at  the  starting  and  ending  depths  are  the 
four  combinations  of  positive  and  negative  angles,  that  is, 

RAY  TYPE  SIGN(0x)  SIGN(0T) 

1  +  + 

2  + 

3  + 

4  - 
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where  0^  is  the  starting  angle  at  the  source  depth,  0y  is  the  ending  angle  at  the  target  depth, 
and  the  signs  are  positive  for  angles  downward  from  the  horizontal.  These  four  ray  types 
are  used  for  calculating  propagation  losses  to  the  target  and  volume  backscatterers.  Ray 
paths  which  terminate  at  turning  points  have  degenerate  ray  types,  so  there  are  only  two 
unique  ray  types  for  paths  to  the  surface  or  bottom: 


RAY  TYPE 

SIGN(0x> 

SIGN(0s) 

S1GN(0b) 

1 

+ 

+ 

- 

2 

_ 

+ 

where  0<j  and  0B  are  the  angles  at  the  surface  and  bottom,  respectively.  These  ray  types 
are  used  for  calculating  surface  and  bottom  reverberation. 

The  ray  tracing  methods  in  the  LORA  computer  program  have  been  revised  to 
reduce  storage  and  to  include  calculations  for  caustics.  Formerly,  in  the  LORA  computer 
program,  rays  were  generated  one  at  a  time,  with  a  stored  ray  path  interpolated  with  the 
current  ray  path  for  target  range.  The  current  ray  path  parameters  were  then  stored  over 
the  old  ray  path  parameters,  a  new  ray  path  for  a  new  starting  angle  was  generated,  etc. 

The  implementation  of  caustic  corrections  necessitated  scanning  families  of 
rays  (rays  with  selected  starting  angles  or  vertex  sound  speeds)  to  isolate  caustics  and 
iterate  on  vertex  sound  speed  to  exactly  determine  caustic  parameters.  With  the  old 
method  of  ray  path  generation,  storage  would  have  been  increased  dramatically. 

In  the  LIRA  program  calculations  are  performed  in  this  order: 

1 .  Ray  segments  are  traced  from  upper  vertexing  depth  to  source  depth,  target 
depth,  and  bottom  depth  foi  40  bottom  bounce  and  40  convergence  zone  (refractive) 
rays.  Five  ray  parameters  are  calculated  for  each  segment:  range,  travel  time,  absorption 
loss,  and  first  and  second  range  derivatives  (with  respect  to  vertex  sound  speed).  The 
total  storage  required  is  2400  cells.  The  storage  areas  required  for  propagation  loss  to  the 
target  and  for  volume  reverberation  are  overlaid  in  core. 

2.  A  backscatterer  (target  is  included  in  this  category)  is  selected. 

3.  A  ray  type  is  selected.  There  are  four  basic  ray  types,  with  rays  either  up  or 
down  at  the  source  or  up  or  down  at  the  target. 

4.  A  ray  family  which  contains  no  discontinuities  in  the  range  derivative  is  selected 
from  the  set  of  40  ray  paths. 

5.  A  cycle  number  (half  the  number  of  ray  reversals)  is  selected.  This  number  often 
corresponds  to  the  number  of  convergence  zones  for  near-surface  source  and  target. 

6.  The  stored  ray  segments  are  used  to  compute  40  ray  paths  for  the  given  cycle 
number  and  ray  type.  The  minimum  and  maximum  range  values  are  inspected  to  ensure 
that  they  overlap  user-specified  ranges  of  interest. 

7.  The  range  derivatives  are  examined  to  find  successive  values  which  straddle 
zero.  Other  tests  are  made  by  comparing  successive  ranges  and  range  derivatives  to 
determine  possible  zeroes  for  the  range  derivative. 

8.  An  iterative  search  is  performed  to  obtain  the  exact  parameters  for  which  the 
range  derivative  is  zero.  Successive  values  for  vertex  sound  speed  Cn  are  estimated  with 
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Cn+1  -  Cn~ 


where  R  is  range. 

9.  The  propagation  loss  is  calculated  for  user-specified  ranges  by  using  Pedersen’s  equations 
(Ref.  2)  adapted  from  Brekhovskikh’s  treatment  (Ref.  4)  for  a  field  near  a  caustic. 
Minimum  and  maximum  ranges  for  the  caustic  field  are  obtained  so  that  ray  interpolation 
in  the  range  region  dominated  by  the  caustic  can  be  excluded. 

10.  If  more  than  one  caustic  occurs  in  this  family,  Steps  8  and  9  are  repeated. 

1 1 .  The  remaining  ray  paths  in  this  family  are  used  to  interpolate  for  propagation  loss  at 
user-specified  ranges  which  lie  outside  all  the  caustic  intervals  of  Step  9. 

12.  The  next  cycle  number  is  selected  and  Steps  6-1 1  are  repeated. 

13.  After  all  the  cycles  have  been  exhausted  that  allow  calculations  for  user-specified  ranges, 
the  next  family  of  rays  is  selected  and  Steps  5-12  are  repeated. 

14.  After  all  the  ray  families  have  been  exhausted,  the  next  ray  type  is  selected  and  Steps 
4-13  are  repeated. 

1  5.  Alter  the  two  or  four  ray  types  have  been  exhausted,  the  next  backscatterer  is  selected 
and  Steps  3-14  are  repeated. 

CAUSTICS 

Detections  at  long  ranges  typically  occur  only  in  the  regions  of  caustics.  Correct 
evaluation  of  caustic  fields  is  therefore  essential  to  estimate  detection  ranges  with  any 
degree  of  confidence. 

A  large  number  of  caustics  can  exist  over  a  range  interval  of  300  nmi.  Usually  only 
one  or  two  main  caustic  peaks  appear  within  the  range  interval  of  a  ray  cycle,  but  often 
enough  there  are  four  to  eight  caustics  per  ray  cycle.  The  deeper  the  source  and  target 
are  below  the  surface,  the  greater  is  the  range  separation  between  caustics  belonging  to 
different  ray  types  in  a  given  ray  cycle  and  the  more  the  caustics  are  resolved  as  individual 
peaks.  A  common  sound-speed  profile  is  illustrated  in  Fig.  1 ,  as  well  as  a  ray  diagram 
(not  to  scale).  A  ray  family  exists  for  which  four  ray  types  with  maximum  sound  speeds 
less  than  the  surface  sound  speed  cross  the  target  depth  at  e,  f,  g,  and  h.  Each  of  the 
four  ray  types  may  produce  a  caustic  at  the  target  depth  (though  not  necessarily). 

Another  ray  family  exists  for  which  four  ray  types  strike  the  surface  and  cross  the  target 
depth  at  a,  b,  c,  and  d.  Again,  each  of  these  ray  types  could  produce  a  caustic  at  the 
target  depth.  For  this  sound-speed  profile,  the  two  ray  families,  each  with  four  ray  types, 
could  produce  eight  caustics  per  ray  cycle.  If  a  ray  cycle  is  60  kyd,  then  a  600-kyd 
plot  of  propagation  loss  could  contain  contributions  from  80  caustics.  A  more  complex 
sound-speed  profile  than  the  one  of  Fig.  I  could  increase  the  maximum  number  of  caustics. 
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Figure  I .  Caustics  possible  at  the  target  depth 
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Figure  2.  The  crossing  of  neighboring  rays. 


For  high  frequencies,  say  above  1  kHz,  standard  ray  tracing  procedures  that  lack 
specialized  treatment  for  evaluating  fields  near  caustics  still  yield  adequate  predictions  for 
detection  ranges.  The  magnitude  of  the  peaks  may  be  in  error,  but  the  ranges  at  which 
signal  excess  crosses  zero  are  reasonably  accurate.  Above  1  kHz  the  high  absorption  of 
sound  limits  active-sonar  detection  to  one  or,  at  most,  two  convergence  zones.  As  frequency 
decreases,  however,  the  peaks  attributed  to  caustic  fields  broaden  and  decrease  in  magnitude. 
Unmodified  ray  theory  cannot  predict  this  broadening  effect,  which  is  nevertheless  known 
from  normal  mode  calculations  for  propagation  loss.  Below  700  Hz  the  peaks  are 
broadened  in  range  by  more  than  a  kyd.  At  500  Hz  absorption  loss  is  only  1  dB  for  500  kyd 
of  path  length.  Below  500  Hz  intermediate-range  propagation  effective  for  detection  appears 
to  be  due  to  the  fields  at  and  near  the  caustics. 


DEFINITION  OF  A  CAUSTIC 


In  ray  theory  a  caustic  occurs  when  neighboring  rays,  i.e.,  rays  contained  within 
some  differential  ray  tube,  coalesce  and  cross  at  some  distance  from  the  source  (see  Fig.  2). 
A  smooth  caustic  is  formed  when  consecutive  neighboring  rays  in  a  finite  ray  tube  cross, 
producing  a  series  of  neighboring  points  constituting  a  caustic  surface  (see  Fig.  3a). 

In  the  geometry  of  the  LIRA  model,  the  vertical  range-depth  plane,  the  smooth  caustic 
becomes  a  line.  Figure  3b  illustrates  schematically  how  two  rays,  finitely  separated,  cross 
near  a  caustic.  As  the  rays  become  arbitrarily  close  together,  point  “p”  moves  onto  the 
caustic  surface,  represented  by  the  dotted  line.  All  rays  forming  the  caustic  are  tangent 
to  it.  Only  the  underside  of  the  caustic  is  illuminated  in  Fig.  3a.  No  rays  of  the  family 
creating  the  caustic  can  pass  into  the  shadow  region.  For  this  reason  the  caustic  is 
sometimes  referred  to  as  an  envelope. 


Ray  theory  predicts  that  a  caustic  intensity  is  zero  in  the  shadow  zone  and  infinite 
on  the  envelope.  Near  the  envelope,  in  the  illuminated  region,  ray  theory  may  predict 
unreasonably  high  intensities,  depending  on  the  interpolation  points. 


In  ray  theory  intensity  per  unit  solid  angle  I/F  is  easily  derived  using  the  concept 
of  solid  angle.  One  of  the  resulting  forms  is 


I/F 


tan  </>x  tan  0Z  CyCz/Cx 


-1 


(1) 


where 

R  is  range, 

C  denotes  sound  speed, 

0  ^enotes  angle  down  from  horizontal, 
and  X,  v,  and  Z  denote  source,  vertex,  and  ending  depths,  respectively. 

The  partial  derivative  is  evaluated  along  the  horizontal,  i.e.,  depth  z  is  held  constant. 
Obviously,  when  rays  cross  at  depth  z,  forming  a  caustic, 


(2) 
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The  intensity  on  the  caustic,  by  Eq.  ( 1 ),  must  be  infinity. 

Actually  the  intensity  on  the  caustic  is  finite.  There  is  also  a  boundary  region  or 
boundary  layer  on  both  sides  of  the  caustic  in  range  and  depth  for  which  ray  theory  is  in 
error.  The  solution  to  the  wave  equation  exactly  describes  the  pressure  field  anywhere  in  the 
space-time  domain.  Unfortunately,  the  wave  equation  has  not,  in  general,  been  solved.  With 
some  restrictions,  the  method  of  normal  modes  produces  exact  solutions  for  special  cases. 

Foi  tunately,  approximation  solutions  have  been  derived  for  the  wave  equation  at  and  near 
caustics. 

APPROXIMATION  FOR  THE  CAUSTIC  FIELD 

Requirements  for  several  scientific  disciplines,  quantum  mechanics,  electromagnetics, 
wave  guide  analysis,  and  acoustics,  have  stimulated  mathematical  development  employing 
normal  mode  theory  (phase  integral  method),  WKB  approximation,  and  uniform  and 
nonuniform  asymptotic  expansion  of  integrals  to  estimate  caustic  fields.  The  solutions 
most  amenable  to  this  work  were  developed  by  Ludwig  (Ref.  3),  whose  expressions  involve 
the  Airy  function  and  its  derivative,  and  by  Brekhovskikh  (Ref.  4),  whose  expressions 
contain  the  Airy  function  alone.  Both  derivations  express  the  field  as  a  superposition  of  an 
infinite  series  of  plane  waves  (as  an  integral),  isolate  the  caustic,  apply  the  method  of 
steepest  descent  (with  different  expressions  for  the  integrand),  and  keep  only  the  significant 
terms  in  frequency  (wave  number).  Brekhovskikh’s  method  has  the  simplification  that  the 
solution  to  the  wave  equation  is  the  WKB  approximation  for  geometrical  optics.  Brekhov- 
skikh's  treatment  was  chosen  for  use  in  the  LIRA  program  because  expressions  in  the  LIRA 
coordinate  system  were  readily  available,  and  the  difference  in  estimations  in  the  region 
of  the  boundary  layer  is  negligible. 

Pedersen's  expressions  (Ref.  2),  derived  from  Brekhovskikh’s  treatment,  give 
propagation  loss  H  as  a  function  of  range  R, 

H(R)  =  -10  log  (1/F)C  -  20  log  |  Ai(t)/Ai(o)] 

(I/F)c  =  2.32  (Cx/Cx  Cv)  cot0xcot0T  f 1 /3 /|R  J  Cy-^§-  2/3> 

/  \  I  3Cv  t 

t  =  (Rc-R)/L 

L  -  0.233!  Cv  rW 

where 

(I/F)c  is  the  intensity  at  the  caustic  relative  to  the  intensity  at  1  yd, 

Ai  is  the  absolute  value  of  the  Airy  function  with  t  as  its  argument, 

Cx  is  the  sound  speed  at  the  source  depth, 

Cy  is  the  sound  speed  at  the  target  depth, 

Cy  is  the  sound  speed  at  the  ray  vertexing  depth, 

0X  is  the  angle  from  the  horizontal  at  the  source  depth. 


(3) 

(4) 

(5) 

(6) 


<t>j  is  the  angle  at  the  target  depth, 
f  is  frequency  in  Hz, 

Rc  is  the  range  of  the  caustic  at  the  target  depth, 


and 


is  the  second  derivative  of  range  with  respect  to  Cy  for  the  ray  path 
starting  at  the  source  depth  and  terminating  at  the  target  depth. 
It  is  evaluated  in  this  case  for  the  value  of  Cv  such  that  =  0 

acv 


For  a  sound-speed  profile  fitted  with  curvilinear  segments,  the  second  derivative 
of  range  with  respect  to  vertexing  sound  speed,  with  depth  held  constant,  is  given  (Ref.  5) 
by 


-  (a  +  5B  +  3R)/Cy 


where 


a  =  2  Dj  (cv  i  j  +  a2i) 
ali  =  B i  j  (cot20j  +  2Dj) 
a2i  =  B2i  (cot20i+i  +  2Dj) 

0R 

B  =  CV9C7=  =[Di<Bli-B2i)-Ril 
R  =  2  Rj 

The  subscript  i  refers  to  the  i^1  layer  and  li  and  2i  refer  to  the  upper  and  lower  boundaries 
of  the  i1'1  layer,  respectively.  Equations  for  Dj,  Bj,  and  B-ij  are  given  in  Ref.  1,  page  22,  or 
in  Ref.  6.  Summations  include  all  segments  needed  to  construct  a  ray  path  with  a  given 
number  of  cycles  and  ray  type  connecting  source  and  target  depths. 

Figure  4  shows  a  plot  of  the  term  -20  log  [Ai(t)/Ai(0)]  used  in  Eq.  (3).  Range  is 
related  to  t  by  Eq.  (5).  The  caustic  envelope,  i.e.,  the  position  at  which  the  range  deri¬ 
vative  is  zero,  corresponds  to  t  =  0.  The  illuminated  side  of  the  caustic  corresponds  to 
negative  t  and  the  shadow  zone  to  positive  t. 

Figure  5  shows  a  plot  of  propagation  loss  near  a  caustic  derived  by  standard  ray 
tracing  and  the  corresponding  propagation  loss  obtained  by  using  Eq.  (3).  The  shape  of 
the  Airy-integral  approximation  is  reversed  from  Fig.  4  to  Fig.  5.  This  is  because  for 
positive  L,  t  is  negative  in  Eq.  (5)  for  R  >  Rc-  The  LIRA  program  uses  Eq.  (3)  in  the 
region 


-2.2  <t  <3 

Outside  this  boundary  layer,  LIRA  uses  standard  ray  tracing  for  derivation  of  propagation 
loss. 

The  LIRA  program  makes  no  attempt  to  smooth  the  changeover  from  Airy-function 

12 


. . . . . — 


t 


Figure  4.  Difference  in  propagation  loss  relative 
to  that  at  the  caustic  envelope  located  at  t  =  0 
(from  Ref.  2). 

caustic  curve  to  the  ray  tracmg  interpolation.  For  small  steps  in  range,  interpolation 
produces  a  small  notch  on  the  illuminated  side  of  the  caustic  (t  =  -2.2),  a  minimum  at 
72.4  kyd  in  the  curve  in  Fig.  6.  For  range  increments  of  1  kyd  or  more,  the  transition 
is  almost  always  smooth. 

The  LIRA  program  does  not  interpolate  the  illuminated  side  of  the  caustic  for 
ranges  at  which  there  are  no  ray  arrivals.  Interpolation  of  a  given  caustic  could  be  truncated 
between  t  =  0  and  t  =  -2.2. 

Locating  a  Caustic 


Locating  the  range  of  a  caustic  for  a  given  target  depth  requires  numerically 
solving  Eq.  (2)  for  Cv  =  Cc  such  that 


D(Cc) 


Since  there  are  many  caustics,  Eq.  (2)  must  be  solved  many  times  to  determine  Cy,  Rc 
and  d^-R/dCy  for  use  in  Eqs.  (4)  and  (6).  It  is  important  not  to  confuse  the  solution  for 
one  caustic  with  that  of  another.  For  a  given  caustic  certain  parameters  are  fixed:  starting 
and  ending  depths,  ray  type,  number  of  ray  cycles,  and  sound-speed  limits  on  ray  family. 

A  sequence  of  ray  segments  is  calculated  for  the  rays  in  the  family  for  a  certain  number 
of  preselected  vertex  sound  speeds  (Cy’s).  The  ray  segments  span  upper  turning  depth 
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Figure  5.  A  caustic  for  unmodified  ray  theory  and  the 
Brekhovskikh  approximation  in  the  boundary  layer. 


Figure  6.  Caustic  interpolation  merged  with  standard 
ray  theory  outside  of  the  boundary  layer  (interpolation 
is  at  0. 1-kyd  intervals). 
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(either  a  refractive  apex  or  surface  hit)  to  source,  target,  and  bottom  depths.  The  ray 
parameters  calculated  for  each  segment  are  range,  range  derivative,  second  range  deriva¬ 
tive,  time,  and  absorption  loss.  The  ray  segments  are  combined  in  the  manner  described 
in  Ref.  1,  pp.  26-30,  to  obtain  the  ray  parameters  for  a  path  of  given  cycle  number  and 
ray  type.  The  range  derivatives  for  the  set  of  rays  in  the  given  ray  family  are  scanned  for 
changes  of  sign.  Whenever  a  sign  change  occurs,  the  range  derivative  must  have  crossed 
zero,  and  a  caustic  has  occurred.  Whenever  a  caustic  is  found,  its  parameters  are  determined 
by  an  iterative  procedure  described  in  the  next  section. 

A  second  scan  of  the  family  of  rays  is  performed  to  try  to  isolate  caustics.  Denote 
iterated  ranges  by  Rj  and  derivatives  by  Dj.  If  Rj+j  <  Rj  and  Dj+|  or  Dj  >  0  or  if  Rj  < 

Rj+1  and  Dj+ j  or  Dj  <  0.  then  two  or  more  caustics  occur  in  the  interval  (Rj,  Rj+j ).  An 
example  is  given  by  Fig.  7.  The  slopes  (Dj,  D->)  are  positive  at  R|  and  R->.  Therefore, 
for  R-i  to  be  less  than  R|,  a  negative  slope  must  occur  in  the  interval  (Rj,  R-,).  The 
derivative  must  pass  through  zero  twice,  First  positive  to  negative,  then  negative  to  positive, 
for  R  to  decrease  between  Ri  and  R->. 


Figure  7.  Range  R  and  range  derivative  D  vs  vertex  sound  speed  Cy-axis. 

In  order  to  isolate  the  two  caustics,  it  is  necessary  to  find  a  point  on  the  Cy-axis 
where  the  derivative  changes  sign.  The  interval  (Cj,  C2)  is  halved  to  obtain  vertex  sound 
speed  C3,  and  another  ray  is  traced  in  order  to  calculate  the  corresponding  range  R3  and 
derivative  D3.  If  the  derivative  has  changed  sign,  then  both  caustics  have  been  isolated, 
one  on  either  side  of  C3,  that  is,  one  each  in  the  intervals  (Cj,  C3)  and  (C3  C2).  If  the 


derivative  has  not  changed  sign.  R  is  compared  to  R  | .  In  Fig.  7  R3  >  R  j ,  so  the  two 
caustics  must  occur  to  the  right  of  C3,the  interval  (C3,  C?)  is  halved,  and  C'3  and  R, 
replace  C,  and  Rj,  respectively.  In  the  LIRA  program,  this  procedure  continues  until  the 
sign  change  in  the  derivative  occurs,  successive  sound  speeds  are  less  than  1 0-5  ft/s  apart, 
or  30  iterations  have  passed.  Typically,  only  one  or  two  iterations  are  needed. 


ITERATION  TO  EXACTLY  DETERMINE  CAUSTIC  PARAMETERS 


Once  a  caustic  has  been  isolated  by  finding  a  sound-speed  interval  (Cj,  C-,]  in  which 
the  range  derivative  changes  sign,  an  iterative  procedure  is  used  to  determine  the  value  C. 
for  vertex  sound  speed  for  which  the  range  derivative  is  zero.  Application  of  the  Newton- 


Raphson  method  gives  successive  approximations  Cn  for  Cc: 


where 


Cn+,  =  Cn-D(Cn)/D'(Cn) 


(7) 


D(Cn)  = 


3R 

dC„ 


and  D'(Cn)  = 

3CJ 


The  estimate  forCc  in  the  interval  [C  j ,  CN]  for  the  first  iteration  is  obtained  by  linearly 
interpolating  between  the  points  (Cj,  Dj)‘and  (CN.  D->)  for  Dc  =  0.  If  Eq.  (7)  yields 
an  estimate  for  Cn  within  f C  j ,  CN],  a  new  ray  is  traced  using  Cn,  and  Eq.  (7)  is  reapplied. 
The  sign  of  D(Cn)  is  compared  to  the  signs  of  D  |  and  D-,  in  order  to  update  the  bracketing 
values  [C ] .  C2]  and  reduce  the  interval  bracketing  the  caustic.  D(Cn)  replaces  D  ]  or  D> 

.  whichever  has  the  same  sign;  correspondingly,  Cn  replaces  Cj  or  C> 

There  are  several  reasons  why  Eq.  (7)  might  fail  to  produce  convergence  to  the 
value  Cc  for  which  D(Cc)  =  0.  The  first  derivative  could  contain  an  extremum  in  the 
interval  [Cj,  C2j  as  in  Fig.  8b.  The  second  derivative  could  be  discontinuous  due  to  the 
discontinuity  of  the  derivative  of  the  sound-speed  gradient  at  the  interface  of  sound-speed- 
profile  segments.  The  caustic  could  be  nearly  horizontal  so  that  the  second  derivative 
approaches  zero.  If  at  any  iteration  Cn  lies  outside  the  bracketing  values  [C] ,  CN],  the 
regula  falsi,  or  secant,  method  is  used.  This  method  is  the  same  as  the  interpolation  to 
obtain  the  initial  estimate  for  Cn-  Figure  8a  shows  the  interpolation  for  a  well-behaved 
rapidly  converging  derivative.  At  each  iteration  the  bracketing  values  [C, , CN]  are  updated 
as  before. 

Sometimes,  Newton’s  method  fails  and  the  regula  falsi  method  converges  too 
slowly.  An  example  of  this  is  shown  in  Fig.  8b.  Successive  iterations  move  Cn  slowly 
to  the  left.  A  maximum  of  30  iterations  is  allowed  per  caustic.  If  the  secant  method  does 
not  produce  a  sign  change  in  the  derivative  or  reduce  the  magnitude  of  the  derivative  by  0.6 
for  successive  derivatives  of  the  same  sign,  then  the  next  Cn  is  selected  to  be  the  midpoint 
of  the  bracketing  interval  [C  j ,  C2).  In  Fig.  8b  halving  the  interval  four  times  will  produce 
a  Cn  for  which  Eq.  (7)  will  again  apply.  On  every  iteration  a  new  attempt  is  made  to  pre¬ 
dict  Cn  using  Eq.  (7).  As  soon  as  Cn  falls  within  the  updated  bracketing  values  [Cj,  C->], 


a.  Rapid  convergence. 


b.  Slow  convergence. 

Figure  8.  Constructions  for  the  regula  falsi  method. 
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the  Newton  method  is  re-employed. 

If  convergence  does  not  occur  within  30  iterations,  the  caustic  parameter  C .  is 
assumed  to  be  the  value  Cn  for  which  the  minimum  absolute  value  of  the  derivative's 
occurred  among  all  the  iterations.  The  corresponding  second  derivative  is  used  to  calculate 
propagation  loss. 
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REVISED  BACKGROUND  NOISE  MODEL 


PROPAGATION  LOSS  FOR  REVERBERATION  CALCULATIONS 

Since  propagation  loss  enters  into  calculation  of  reverberation,  the  changes  in  the 
propagation  loss  model  described  in  the  previous  section  impact  the  calculation  of  reverbera¬ 
tion.  The  propagation  loss  model  used  for  reverberation  calculation  is  the  same  as  that  used 
for  propagation  loss  to  the  target.  Target  depth  is  replaced  by  surface,  volume,  and  bottom 
depths  for  surface,  volume,  and  bottom  reverberation  calculations,  respectively.  The  depth 
used  for  the  volume  backscatterers  is  one-half  of  the  user-specified  column  depth.  This 
effectively  converts  the  scattering  column  to  a  scattering  layer.  The  set  of  40  bottom 
bounce  and  convergence  zone  (refractive)  rays  is  different  from  that  used  in  calculating 
propagation  loss  to  the  target,  because  target  depth  enters  into  the  selection  of  the  rays. 

It  is  possible,  if  one  set  of  rays  were  used  for  both  cases,  to  change  reverberation  results 
by  changing  the  target  depth.  Therefore,  target  depth  is  replaced  by  volume  reverberation 
depth  and  a  new  set  of  80  rays  is  used  for  reverberation  calculations. 

PULSE-AVERAGED  REVERBERATION 

The  LIRA  computer  program  allows  a  maximum  pulse  length  of  100  s.  For  this 
discussion  a  pulse  length  t  that  is  1  <  r  <  100  is  considered  to  be  a  long  pulse  length. 

The  LORA  model  allows  only  short  pulse  lengths.  In  LORA  the  form  of  the 
reverberation  is  (Ref.  1 ) 

LK  =  Lq  +  2Hk  +  Sk  +  10  log  (Ak/Ad)  (8) 


where 

L()  is  the  source  strength. 

H|^  is  the  propagation  loss  including  beam  pattern  losses, 

Sj(  is  the  backscattering  strength  per  AQ, 

Aj^  is  the  backscattering  area, 
and  AQ  is  the  reference  area. 

For  short  pulse  lengths  Hj^  does  not  vary  over  the  area  A^.  For  long  pulse  lengths,  however, 
may  vary  considerably  over  the  range  dimension  of  ak.  In  the  case  of  an  FM  or  PRN 
signal  and  matched  filter  processing,  range  is  highly  resolved.  For  a  CW  pulse  length  of  r,  the 
range  interval  used  to  compute  reverberation  area  is  Cr/2  cos  0<j.  For  long  pulse  lengths 
and  a  CW  signal,  propagation  loss  for  reverberation  is  averaged  over  this  range  interval. 

After  matched  filter  processing,  the  area  with  FM  or  PRN  signals  is  Ct/2  Wr  cos  = 
C/2  W  cos  For  large  rW  products,  the  range  resolution  is  small.  For  a  10-s  pulse  and 
100-Hz  bandwidth,  for  example,  the  range  interval  contributing  to  Aj^  is  25  ft.  is 
assumed  not  to  vary  over  this  interval,  and  Eq.  (8)  applies  to  the  LIRA  signal  waveforms. 

In  the  LIRA  model  this  reverberation  is  referred  to  as  “instantaneous  reverberation.” 

An  optional  plot  of  “Instantaneous  Reverberation  vs  Target  Range”  is  one  of  the  program 
outputs. 
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The  instantaneous  reverberation  corresponds  to  the  time  when  the  leading  edge  of 
the  target  echo  arrives  at  the  receiver.  Since  the  signal  is  correlated  for  the  duration  of  the 
pulse  length,  instantaneous  reverberation  must  be  averaged  over  that  interval;  that  is, 


L(T2,  r) 


1  T  ~>  +  T 

10  log  [Lf  ~ 

7  '2 


1001  L<l)dt] 


(9) 


where  T2  is  two-way  travel  time  to  the  target,  and  L  is  reverberation  level  in  dB.  Thus, 
averaged  reverberation  is  obtained  by  sliding  an  averaging  window  of  length  r  over  the 
function  of  instantaneous  reverberation  vs  time.  This  averaging  is  performed  in  the  LIRA 
program  for  surface,  bottom,  and  volume  reverberation  for  each  two-way  time  to  the  target 
corresponding  to  a  user-specified  target  range.  The  total  pulse-averaged  reverberation  is 
the  incoherent  sum  of  the  bottom,  surface,  and  volume  reverberation  components.  The 
plot  “Averaged  Reverberation  vs  Range”  is  an  optional  program  output. 

Reverberation  is  computed  only  for  those  target-range  values  (and  corresponding 
time  values)  entered  by  the  user.  The  reverberation  pressure  squared  is  considered  to  lie  on 
straight  lines  connecting  these  range  (time)  values.  Therefore,  Eq.  (9)  is  the  integral  of 
reverberation  as  a  piecewise,  linear  function  of  time.  If  the  user  neglects  to  enter  a  closely 
spaced  sequence  of  range  values,  the  piecewise,  linear  reverberation  function  could  be 
greatly  in  error.  Therefore,  for  values  of  pulse  length  which  will  require  integration,  that 
is  for  r  >  1  s,  the  user  is  encouraged  to  enter  his  range  array  incremented  in  steps  of  2  kyd 
or  less.  Since  instantaneous  reverberation  values  are  computed  only  for  range  values  that 
the  user  enters,  the  pulse  averaging  near  and  at  the  end  of  the  sequence  of  ranges  will  not 
extend  over  a  whole  pulse  length.  In  particular,  the  last  value  in  the  pulse-averaged  rever¬ 
beration  array  will  be  the  value  of  the  instantaneous  reverberation  at  the  target  range.  The 
averaging  at  the  end  is  ^ 

UT2,Te)=  lOlog^  f  10ai  L(t)dtj.r>TE-T2 

where  TE  is  the  time  corresponding  to  the  ending  range.  This  “edge  effect”  in  averaging 
appears  between  target  ranges  RE  -0.83t  and  Rp  kyd,  where  Rp  is  the  largest  range  speci¬ 
fied  by  the  user.  Thus,  a  10-s  pulse  would  have  incorrect  averaging  for  the  last  8.3  kyd 
(1/2  range  covered  by  t  =  10  s)  before  the  maximum  range,  and  a  100-s  pulse,  83  kyd. 

The  user  must  be  sure  that  his  input  maximum  range  exceeds  his  maximum  range  of  interest 
by  0.83t.  This  correction  applies  if  a  pulse  length  greater  than  1  s  is  used  and  if  the  user  is 
concerned  about  time-averaged  vs  instantaneous  reverberation  values. 
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SEPARATION  OF  SOURCE  AND  RECEIVER  IN  DEPTH 


The  LIRA  program  allows  the  source  and  receiver  to  be  separated  in  depth  (one 
dimension)  but  not  in  range  or  azimuth  (two  or  three  dimensions).  Propagation  loss, 
travel  time,  and  reverberation  for  a  given  target  range  are  obtained  by  averaging  these 
quantities  computed  for  source/receiver  at  the  source  depth  and  source/receiver  at  the 
receiver  depth.  Justification  of  this  procedure  and  estimates  of  the  error  introduced  into 
the  reverberation  calculations  constitute  the  remainder  of  this  section.  The  method  em¬ 
ployed  in  LIRA  is  convenient  because 

•  The  complex  computational  and  data  storage  structure  is  essentially  the  same 
as  for  the  case  of  source  and  receiver  at  the  same  depth. 

•  Computation  time  is  drastically  reduced  for  what  is  otherwise  a  degenerate  case 
of  a  bistatic  geometry. 

•  Data  storage  is  correspondingly  reduced. 

PROPAGATION  LOSS 

The  propagation  loss  for  the  principal  ray  (largest  arrival)  from  the  source  to  the 
target  to  the  receiver  is 

H  j  2  =  H,  +  H2  =  (2H,  +2H2)/2 

where  the  subscripts  1  and  2  refer  to  the  source-target  and  target-receiver  paths,  respectively. 
The  travel  time  for  the  source-target-receiver  path  is 

T12  =  T1  +T2  =  <2T,  +2T2)/2 


Therefore,  propagation  loss  for  the  source-target-receiver  path  can  be  computed  by  cal¬ 
culating  two-way  propagation  loss  for  the  source-receiver  at  the  source  depth,  again  at  the 
receiver  depth,  and  averaging.  The  same  can  be  done  with  time.  No  inherent  error  is 
introduced  into  propagation  loss  or  time  by  using  this  procedure. 

The  same  arguments  apply  to  propagation  loss  derived  from  the  incoherent  sum  of 
the  multipath  intensities.  Propagation  loss  from  the  source  to  the  target  is 

HX  =  -10  log  (E  10"0-1  HL) 
i 

and  from  the  receiver  to  the  target  is 

HR  =  -10  log  (2  1 0_0  1  HJ) 

j 

The  propagation  loss  for  a  source-target-receiver  path  is 

H^-IOIogllO-OKHW), 
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The  total  propagation  loss  from  all  source-target-receiver  paths  is 


which  is  equivalent  to 

H  =  -10  log  £  1(T01H*  -10  logS  10-01  Hlj 
i  J 

*  HX  +HR 
=  (2HX  +  2HR)/2 

Therefore,  the  total  propagation  loss  for  all  of  the  incoherently  summed  source- target- 
receiver  paths  is  equivalent  to  the  average  of  the  two-way  propagation  losses  from  source 
to  target  and  from  receiver  to  target. 

REVERBERATION 

In  LIRA  the  reverberation  level  is  evaluated  for  the  two-way  time  to  target,  assum¬ 
ing  the  path  is  source-target-source.  Again,  reverberation  level  is  computed  for  the  receiver- 
target-receiver  path  as  if  the  receiver  were  a  source.  The  reverberation  levels  for  a  given 
target  range  are  averaged,  and  the  result  is  an  approximation  for  reverberation  for  separated 
source  and  receiver.  The  following  discourse  is  a  justification  for  using  the  two-pass  approxi¬ 
mation  for  reverberation  when  source  and  target  are  separated  in  depth. 

TRUE  REVERBERATION 

Quantitatively,  reverberation  level  for  a  given  target-echo-return  time  for  source 
and  receiver  separated  in  depth  is  expressed  by 

L  =  I  Lk 
k  K 

Lk  =  Lo  +  BX^X*  +  BR^R^  "  Hk  “  Hk  +  S]^^,  0^)  +  10  log  (Aj*R/AQ)  (10) 

where 

LQ  is  the  source  level, 

k  is  the  index  of  the  backscattering  element, 

B^  is  the  source  beam-pattern  response  for 

<t>\  the  ray  angle  with  respect  to  the  horizontal  of  the  source, 

Br  is  the  receiver  beam  pattern  response  for 
0r  the  ray  angle  at  the  receiver, 

H|*  is  the  one-way  propagation  loss  between  source  and  the  k*’1  backseat terer, 

HR  is  the  one-way  propagation  loss  between  the  receiver  and  the  k1*1  backscatterer. 


Sfc  is  the  bistatic  backscattering  strength  in  dB  per  unit  area,  where  ^  and  <t>l  are 
angles  of  incidence  and  scattering  (return),  respectively, 

Aq  is  a  reference  area  which  has  the  same  spatial  units  as  the  backscattering  area, 
and 

X  R  th 

Ak  is  the  effective  scattering  area  of  the  ktn  scattering  surface  or  cross-sectional 

areas  associated  with  column  volume  scattering. 

TWO-PASS  APPROXIMATION  FOR  REVERBERATION 

The  two-pass  method  computes  reverberation  Lj*  for  the  source  at  source  depth, 
reverberation  for  the  source  at  receiver  depth,  and  averages  the  result. 

Lk  =  Lo  +  BX^X^  +  BX^X^  “  2Hi?  +  ^k^'k'  4>ic)+  10  lo8  (Ak  /A0) 

Lk  =  Lo  +  BR^R^  +  BR^R)_  2Hk  +  Sk^k’  $[)  +  10  loS  (Ak  /A0) 

Adding  and  dividing  by  2  gives  the  reverberation  level  Lk  R  for  the  kth  backscatter 
obtained  by  the  two-pass  method: 

Lk  R  =  +  BX(0X)  +  BR(0R)  -  Hf  -  H* 

+  1/2  [Sk(4,  0>k)  +  Sk(4,  0[)]  +  10  log  [(A^A^)>/2/A0]  (11) 

The  true  reverberation  level  for  separated  source  and  receiver  is  given  by  Eq.  (10):  the 
reverberation  predicted  by  LIRA  is  given  by  Eq.  (11).  The  two  equations  differ  only  in 
the  scattering  strengths  and  backscattering  areas. 

ERROR  IN  THE  BACKSCATTERING  AREA  FOR  THE  TWO-PASS  METHOD 

True  Backscattering  Area.  The  true  backscattering  area  is  illustrated  in  Fig.  9.  Paths 
a  and  b  are  incident  upon  the  surface  with  angle  6r  The  difference  in  path  length  between 
paths  a  and  t  is  C  (r|  +  where  r  =  r  j  +  r2  is  the  ping  length.  All  points  on  the 
surtace  over  length  Ar  will  make  a  scattering  contribution  to  the  reverberation  at  the 
receiver  at  time  T2,  the  time  at  which  the  leading  edge  of  the  target-echo  pulse  arrives  at 
the  receiver.  Path  b  is  traveled  by  the  leading  edge  of  the  transmitted  pulse,  path  a 
the  trailing  edge. 

The  backscattering  area  for  the  geometry  of  Fig.  9  is  r  A  0  A  r,  where  r  is  the  hori¬ 
zontal  range  from  source/receiver  to  the  backscattering  area,  and  A  0  is  the  azimuthal 
angle  containing  significant  source  energy.  The  true  backscattering  area  is  derived  by 
calculating  the  value  for  Ar.  From  Fig.  9 

Ar  =  Crj  sec  =  Cr-tsec0r 

where 


r  =  r  |  +r2. 


23 


Therefore, 

Cr  sec  0 j  sec  0r 

A  r  = - - - — 

sec  0 j  +  sec  0r 

The  true  backscattering  area  to  be  used  in  Eq.  (10)  is  therefore 

A***  =  rA0Ar  =  r  A  0  C  r  sec  0j  sec  0r /(sec  0j  +  sec  0r)  (12) 

Two-pass  Backscattering  Area.  With  source  and  receiver  separated  in  depth,  LIRA 
calculates  reverberation  as  if  the  source  and  receiver  were  both  located  at  source  depth  and 
then  as  if  both  source  and  receiver  were  at  receiver  depth  and  averages  the  results.  The 
geometry  of  Fig.  9  does  not  apply  to  this  case.  Instead  the  two  passes  require  the  geometry 
of  Fig.  10  to  explain  the  reverberation  backscattering  area. 

The  two  areas  of  Eq.  ( 1  1 )  are 

A*  =  rA0Arj  and  Aj^  =  rA0A^ 

Since  Arj  =  1/2  Cr  sec  0j  and  Ars  =  1/2  Cr  sec  0r  the  geometric  mean  used  in  Eq.  ( 1 1 ) 
is 

(Aj*  Aj^)’/^  =  r  A  0  (Ar j  A^)*^  =  l/2r  A  0  C  r  (sec  0j  sec  0rW“ 
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Figure  10.  Monostatic  backscattering  geometry  for:  (a)  ray  path  from  source  depth;  (b)  ray  path  from 
receiver  depth. 


This  approximated  area  is  to  be  compared  with  the  true  reverberation  area  of  Eq. 
(12)  used  in  Eq.  (10).  The  error  in  reverberation  level  EA  due  to  approximating  back- 
scattering  area  with  the  two-pass  method  is 


Obviously,  when  source  and  receiver  are  at  the  same  depth,  =  <pr  and  EA  =  0.  For  shallow 
angles  of  incidence  and  return,  EA  is  also  very  nearly  zero.  The  greatest  discrepancy  occurs 
when  one  angle  is  small  and  the  other  is  large.  For  an  extreme  case,  say  0j  =  0  and  0r  = 

60  deg,  Ea  =  0.26  dB.  Only  for  very  short  ranges,  in  the  order  of  1  kyd,  could  the  transmit 
and  return  paths  of  significance  be  greatly  different;  yet  even  if  they  are,  the  error  in  approxi¬ 
mating  backscattering  area  is  small  when  using  the  two-pass  method. 

ERROR  IN  BACKSCATTERING  STRENGTH  FOR  THE  TWO-PASS  METHOD 

The  true  backscattering  strength  used  in  Eq.  (10),  Sk  (</>'k,  </>£),  is  bistatic.  The 
approximated  backscattering  strength  used  in  Eq.  (1 1)  for  the  two-pass  method 

2  fSk  <4-  +  \  1 

is  the  average  of  two  monostatic  backscattering  strengths,  i.e.,  the  angle  of  incidence  equals 
the  angle  of  return.  The  error  in  backscattering  strength  Eg  in  dB  introduced  in  reverbera¬ 
tion  level  by  using  the  two-pass  method  is 


Es  =  \  [sk  (4  4) +  sk  (^ic>  4)  ]  -  sk  (4-  4)  0 1 3) 

Obviously,  for  any  backscattering  strength  which  is  isotropic,  Eg  =  0.  Therefore,  volume 
reverberation  predictions  suffer  no  error  from  the  two-pass  method. 

Error  in  Surface  Backscattering  Strength.  Bistatic  surface  backscattering  is  much 
less  understood  than  monostatic.  Most  experiments  interpret  backscattering  strengths 
monostaticaliy.  This  author  knows  of  only  one  functional  form  for  bistatic  surface  back- 
scattering  strength  [after  Brown  and  Saenger,  (Ref.  7)]: 

S(0j,  0r)  =  Aj  +  A2  (sin  +  sin  0r)2  (14) 

where  0j  and  </>r  are  angles  of  incidence  and  scatter,  respectively,  and  Aj  and  A->  are  con¬ 
stants  derived  from  data.  For  this  analysis  only  A->  is  needed.  From  the  data  of  Brown 
and  Saenger,  which  spans  frequencies  from  62.5  to  5000  Hz  and  angles  20  deg  <0j  <  70  deg 
and  1.5  deg  0r  <  70  deg,  A2  has  a  maximum  value  of  10  dB  for  a  weakly  shadowed  sur¬ 
face  and  1 5  dB  for  a  strongly  shadowed  surface. 

The  error  Eg  in  evaluating  the  backscattering  strength  using  the  two-pass  (monostatic) 
method  is  derived  by  substituting  Eq.  (14)  into  Eq.  (13). 

ES  =  4  1^2  (2  sin  0j)“  +  A2  (2  sin  <pT)~]  -  A2  (sin  0j  +  sin  0r)“ 

=  A2  (sin  sin  <pT)~  (15) 

Figure  1 1  shows  the  maximum  error  for  maximum  angle  differences,  i.e.,  one 
angle  is  0  deg  and  the  other  is  at  the  maximum.  Figure  1 2  shows  the  error  if  0j  and  0r 
are  not  more  than  10  deg  apart.  These  two  figures  are  based  on  Eq.  (15).  Errors  for  values 
of  0j  <  20  deg  are  not  computed  because  of  the  limits  of  the  data  upon  which  the  bistatic 
function  is  based.  For  ranges  at  one  convergence  zone  and  greater,  angles  of  incidence  and 
return  are  limited  by  the  bottom;  they  are  less  than  15  deg  and  tend  to  be  about  5  deg. 

Figure  1 1  indicates  that  for  very  short  ranges  at  steep  angles,  significant  error  could  be 
induced  in  surface  backscattering  strength  if  the  two-pass  method  is  used.  For  long  ranges, 
however,  the  worse  case  would  be  less  than  1.8  dB.  Figure  12  shows  that  if  the  angles 
are  not  greater  than  10  deg  apart,  the  maximum  error  is  0.42  dB.  It  is  very  unlikely  that 
this  spacing  in  bistatic  angles  would  be  exceeded  for  any  source/receiver  configuration  in 
any  ocean.  Much  more  likely,  the  angles  will  be  within  2  deg  of  each  other  with  negligible 
induced  error. 

Error  in  Bottom  Backscattering  Strength.  The  Lambert’s  law  approximation  for 
bottom  backscattering  strength  is  of  the  form  [after  MacKenzie  (Ref.  8)]; 

Sr  =  B  |  +  1 0  log  (sin  0j  sin  0r) 

where  0 j  and  0r  are  angles  of  incidence  and  return,  respectively,  and  B  j  is  a  constant.  The 
error  Eg  induced  in  reverberation  level  by  using  the  two-pass  method  in  place  of  the  bistatic 
geometry  is 
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EB  =  T  f 10  l°S  (sin2  <t> i)+  10  log  (sin2  0r)J  -  10  log  (sin  0;  sin  0r) 

=  0 

Therefore,  where  Lambert  s  law  applies,  the  two-pass  method  introduces  no  error  compared 
with  the  bistatic  model.  For  the  steeper  angles,  where  specular  facet  reflection  dominates, 
Lambert’s  law  does  not  hold.  However,  angles  of  incidence  and  return  will  be  nearly 
equal  even  though  source  and  receiver  are  quite  separated  in  depth  and  the  bistatic 
geometry  approaches  monostatic. 
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RESULTS 


In  Appendix  D  LIRA  is  compared  with  three  other  models,  PLRAV,  FACT  and 
AP2,  with  regard  to  propagation  loss.  For  most  cases  all  the  models  agree  well.  Differences 
are  usually  explainable.  The  importance  of  the  comparison  is  to  show  that  LIRA  produces 
a  reasonable  propagation  loss  curve.  LIRA’S  propagation  loss  examples  exhibit  no  path¬ 
ological  excursions  from  the  norm  ct  the  other  models,  nor  trends,  nor  biases.  LIRA'S 
caustic  correction  algorithm  performed  without  idiosyncracies. 

Reverberation  is  calculated  by  means  of  the  same  algorithm  used  for  propagation  loss. 
Because  of  the  way  in  which  reverberation  is  represented,  it  is  difficult  to  make  meaningful 
comparisons  with  other  models  and  with  experimental  measurements. 

The  LIRA  program  has  been  executed  thousands  of  times  and  is  a  reliable  per¬ 
former.  Inputs  and  outputs  are  user  oriented,  designed  for  clarity.  Active-sonar  performance 
predictions  for  targets  from  1  to  600  kyd  from  the  source  cost  $2  to  SI  5  using  a  UNI  VAC 
1 1  JO  computer,  with  plots  costing  extra  (about  $4/plot).  The  four  example  runs  in  Appen¬ 
dix  E  cost  $25,  with  10  min  ot  cpu  time  used.  The  14  propagation  loss  curves  for  the 
cases  in  Appendix  D  cost  S4.70  and  used  107  s  of  cpu  time.  Thus,  one  case  averaged 
$0.34  for  7.64  s  of  cpu  time. 

LIRA  uses  57,000  words  of  computer  core  storage.  The  program  must  be  seg¬ 
mented  (overlaid  in  core)  in  order  to  incorporate  the  plotting  routines. 
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SUMMARY  OF  MODEL  AND  SOFTWARE  IMPROVEMENTS 

LIRA  surpasses  LORA’s  capabilities  in  the  following  ways: 

1 .  Allowable  frequencies  are  from  25  Hz  to  25  kHz. 

2.  Maximum  range  is  1000  kyd. 

3.  Brekhovskikh’s  caustic  corrections  are  applied  to  propagation  losses  to  the 
target  and  all  backscattering  elements  which  produce  reverberation. 

4.  Linear  FM  and  PRN  signals  are  modeled  by  averaging  the  reverberation  time 
sequence  over  the  pulse  length. 

5.  The  source  and  receiver  are  not  restricted  to  being  at  the  same  depth. 

Outputs  are  modified  in  the  following  ways: 

1 .  Calcomp  plots  can  be  made  for  propagation  loss,  reverberation  (including 
components  from  bottom,  surface,  and  volume  backscattering),  and  signal  excess.  These 
are  displayed  vs  range  from  source  to  target. 

2.  Target  ranges  are  user  specified.  Formerly,  the  user  specified  the  number  of 
convergence  zones  and  bottom  bounces,  and  the  output  included  all  ranges  encompassed 
by  the  ray  paths  generated. 

3.  The  output  is  tabulated  for  the  user-specified  target  ranges.  Formerly  the 
output  heading  described  the  mode  of  propagation  (i.e.,  direct  path,  convergence  zone, 
bottom  bounce,  etc.),  and  the  output  was  broken  up  into  various  zones. 

New  inputs  available  to  the  user  are: 

1.  Arrays  of  surface  and  bottom  backscattering  strengths  vs  grazing  angle. 
(Formerly,  an  internal  model  had  to  be  used.) 

2.  Array  of  target  ranges. 

3.  A  new  array  of  flags  for  special  intermediate  printout  for  debugging  purposes. 

4.  Array  of  recognition  differential  vs  reverberation-to-noise  ratio. 

5.  Flag  to  turn  on  or  off  the  caustic  corrections. 

6.  Flag  to  suppress  any  or  all  portions  of  printing  or  plotting. 

7.  Flags  to  exclude  up-  or  down-going  rays  at  the  source  and/or  the  receiver 
(to  model  the  array  on  a  steeply  sloping  bottom). 

8.  Flags  to  select  CW,  FM,  or  PRN  signal  waveform. 

9.  Receiver  depth  (when  different  from  source  depth). 
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RECOMMENDATIONS 

LIRA  experiences  the  same  difficulties  with  submerged  ducts  that  LORA  did. 
Study  should  be  undertaken  to  develop  or  choose  a  model  which  corresponds  in  normal 
mode  theory  to  weak  trapping  of  modes. 

A  big  step  toward  true  bistatic  simulation  of  active  sonar  is  to  allow  source  and 
receiver  to  be  separated  in  the  range  dimension.  With  this  modification  LIRA  would  have 
source  and  receiver  separated  in  two  of  the  three  dimensions  (azimuth,  the  remaining 
dimension). 

Representation  of  reverberation  at  the  output  of  the  signal  processor  causes 
awkwardness  in  comparing  LIRA’S  reverberation  with  that  of  other  reverberation  models. 
A  plot  of  normalized  reverberation  vs  time  would  be  useful. 

Fields  near  horizontal  caustics  cannot  be  represented  by  Brekhovskikh’s  methods. 
A  study  should  be  performed  to  derive  the  fields  for  these  cases.  The  equations  for  fields 
near  horizontal  cusped  caustics  are  available  (implemented  in  the  FACT  model)  (Ref.  9). 
These  equations  could  be  included  in  LIRA,  but  the  rather  rare  occurrence  of  horizontal 
cusped  caustics  (source  and  target  must  be  at  the  same  depth)  reduces  the  priority  of 
this  effort. 
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LOGICAL  STRUCTURE  OF  THE  LIRA  COMPUTER  PROGRAM 

The  LIRA  computer  program  consists  of  two  segments,  one  for  calculations  and  one 
for  plotting.  The  MAIN  routine  (See  Fig.  A-l)  calls  SUBROUTINE  LIRA,  which  causes 
segment  1  to  be  loaded.  It  is  in  this  segment  that  the  user  enters  his  inputs  in  SUBROUTINE 
READIN.  After  all  of  the  user-specified  calculations  have  been  completed  and  if  plot  flags 
have  been  set,  control  returns  to  MAIN,  and  segment  2  is  loaded.  All  storage  areas  except 
those  common  to  both  segments  are  lost.  The  appropriate  plots  are  created  in  segment  2, 
and  control  is  returned  to  the  user  in  READIN  after  passing  through  MAIN  to  reload 
segment  1. 

All  inputs  are  entered  through  SUBROUTINE  READIN.  The  input  parameters  are 
flagged  to  determine  if  they  belong  to  sound-speed-profile  (SSP),  environmental,  reverberation, 
or  detection  categories.  Depending  on  the  inputs,  redundant  calculations  are  skipped  on 
successive  runs  performed  during  the  same  program  execution. 

If  a  user  on  a  demand  terminal  executes  the  program,  and  after  a  given  run  decides 
that  he  would  like  to  plot  his  previous  calculations,  he  can  set  the  appropriate  plot  flags, 
type  the  data  literal  “GOPL,”  and  the  segment  swapping  and  plotting  will  be  performed. 
Control  returns  to  the  user  in  READIN  in  segment  1,  with  the  input  values  the  same  as 
before  he  typed  “GOPL.”  He  can  continue  to  make  runs  with  his  data,  but  unless  he  changes 
the  plot  flags,  plotting  will  still  be  performed  after  each  run. 

The  flag  ISPLIT  is  initialized  as  0.  If  source  depth  ZX  and  receiver  depth  ZR  are 
different,  ISPLIT  is  set  to  1  for  the  first  pass  for  calculating  propagation  loss  and  reverbera¬ 
tion.  For  the  first  pass  ZR  replaces  ZX  (ZX  is  temporarily  stored  as  ZXX).  The  first-pass 
propagation  loss  and  reverberation  arrays  are  stored.  ISPLIT  is  set  to  2.  The  second  pass 
is  performed  with  ZX  restored  to  the  original  source  depth  (ZXX).  Propagation  loss  and  rever¬ 
beration  arrays  are  combined  after  both  passes.  ISPLIT  is  reset  to  0. 

Special  outputs  are  allowed  for  propagation  loss  and  reverberation  if  the  user  enters 
values  for  the  flags  OU TP  and  OUTR,  respectively.  If  OUTP  =  1 ,  all  reverberation  and  detec¬ 
tion  calculations  are  skipped  and  two-way  travel  time  to  the  target,  angle  of  the  largest 
arrival,  the  propagation  loss  of  the  largest  arrival,  and  the  propagation  loss  for  the  incoherent 
sum  of  all  of  the  ray  arrivals  are  printed  for  all  of  the  user-specified  target  ranges. 

If  OUTR  =  1  or  2,  bottom,  surface,  and  volume  components  of  reverberation  as  well 
as  contributions  from  each  previous  ping  are  printed  for  the  times  corresponding  to  two- 
way  travel  time  for  the  user-specified  target  ranges.  If  OUTR  =  1,  the  detection  calculations 
are  performed;  if  OLfTR  =  2,  they  are  skipped. 

Unless  the  user  enters  NOPR,  the  LIRA  program  prints  detection  outputs  under 
the  following  headings: 

•  RANGE,  all  user-specified  target  ranges  in  kyd. 

•  TIME,  two-way  travel  time  for  the  ray  which  gives  minimum  loss  to  the  target. 

•  LOSS,  minimum  two-way  propagation  loss  (or  incoherent  sum),  including  beam 
pattern  corrections. 


33 


•  ANGLE,  angles  at  source  depth  for  rays  which  give  minimum  loss  to  the 
target. 

•  SIGNAL,  level  of  the  return  signal  at  the  output  of  the  signal  processor. 

•  REVERB,  reverberation  level  at  the  output  of  the  beamformer  (corrected  for 
pulse  averaging)  or  at  the  output  of  the  signal  processor  (corrected  for  target 
doppler,  pulse  averaging,  and  processing  gain),  depending  on  whether  recogni¬ 
tion  differential  or  detection  threshold  is  used. 

•  REV  LIM,  “YES”  to  indicate  if  reverberation  exceeds  noise  by  more  than 
3  dB;  “NO”  otherwise. 

•  S/MN,  signal-to-masking-noise  ratio  in  dB. 

•  EXCESS,  signal  excess  in  dB,  with  zero  corresponding  to  50%  probability 
of  detection  if  recognition  differential  is  used  or  the  probability  of  detection 
for  a  given  detection  threshold  if  that  is  used. 

•  PROB  DETECT,  probability  of  detection,  with  0.5  corresponding  to  signal 
excess  =  0  (not  applicable  for  other  probabilities  of  detection  at  signal  excess 
=  0). 
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Figure  A-l  (a).  Flow  chart  for  LIRA  input,  output,  and  calculations. 
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Figure  A- 1  (h).  Flow  chart  for  l.IRA  input,  output  and  calculations. 
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INPUTS  TO  THE  LIRA  COMPUTER  PROGRAM 

Table  B-l  contains  a  list  of  the  most  commonly  used  inputs  to  the  LIRA  computer 
program.  Following  Table  B-2,  the  inputs  are  described  in  greater  detail  along  with  other 
inputs  not  needed  by  the  general  user. 

The  column  under  “Name”  in  Table  B-l  contains  the  name  of  the  variable  or  array 
to  be  typed  by  the  user  when  he  enters  data  according  to  the  format  described  in  Appendix 
C.  Some  variables  have  optional  names.  The  column  under  “Description”  contains  brief 
definitions  of  the  input  parameters.  “Range”  is  the  range  of  values  a  given  parameter  is 
allowed.  Any  value  outside  this  range  is  either  set  to  a  “Typical  Value”  or  calculated  in 
the  program. 

Table  B-2  contains  a  list  of  data  literals  used  to  exert  specific  control  over  the 
program. 

Not  all  of  the  input  parameters  need  to  be  entered  by  the  user.  The  only  input 
required  tor  the  program  to  run  is  sound-speed  profile  (SSP).  Typical  or  calculated  values  are 
used  for  every  input  parameter  not  entered  by  the  user. 

Input  parameters  fall  into  three  categories,  environmental,  reverberation,  and  detection. 
The  computational  structure  of  the  LIRA  program  may  be  thought  of  as  three  nested  loops, 
with  the  environmental  loop  on  the  outside  and  the  detection  loop  on  the  inside.  During  a 
single  program  execution,  many  runs  can  be  made.  On  the  first  run,  unless  otherwise  speci¬ 
fied  by  the  user  through  the  inputs,  the  environmental,  reverberation,  and  detection  loops 
are  all  performed.  Thereafter,  only  those  loops  are  performed  that  are  required  by  category 
of  inputs  entered  by  the  user.  If  he  enters  only  detection  parameters,  only  the  detection 
loop  is  performed.  If  he  enters  reverberation  parameters,  the  reverberation  and  detection 
loops  are  performed.  Savings  in  computer  time  and  cost  are  obtained  by  looping  on  detection 
and/or  reverberation  parameters  rather  than  environmental  parameters.  Table  B-l  groups 
the  input  parameters  into  environmental,  reverberation,  and  detection  categories  so  that 
the  user  can  choose  the  order  to  enter  his  data  and  reduce  his  computer  costs. 

ADDITIONAL  DESCRIPTION  OF  INPUTS 

Certain  inputs  need  to  have  their  descriptions  elaborated  to  clarify  their  meanings 
or  to  acquaint  the  user  with  additional  options: 

BD,  ZBM.  If  bottom  depth  (ft)  is  not  entered,  it  is  assumed  to  be  the  maximum 
depth  in  the  sound-speed  profile.  If  BD  is  entered,  the  sound-speed  profile  is  either  linearly 
interpolated  or  extrapolated  to  find  the  sound  speed  for  the  bottom  depth. 

FREQ,  F.  Frequency  may  be  entered  in  units  of  Hz  or  kHz.  The  program  tests  the 
value  to  determine  the  units.  If  25  <  FREQ  <  100,000,  the  program  assumes  Hz.  If 
0.0001  <  FREQ  <  25,  the  program  assumes  kHz. 

R,  RANGE.  Range  of  target  in  kiloyards;  an  array  of  ranges  and  range  increments 
for  which  propagation  loss,  reverberation,  signal  excess,  and  probability  of  detection  are 
calculated.  The  elements  of  array  R  are  scanned  in  order  to  generate  the  range  array  internal 
to  the  program.  Each  successive  set  of  three  elements  in  R  is  examined  to  see  if  the  third 
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Table  B-l.  Most  commonly  used  LIRA  inputs. 
ENVIRONMENTAL  PARAMETERS 


Name 

Description 

Units 

Range 

Typical  Value;  Comments 

SSP 

Sound-speed  profile.  Array 

ft,  ft/s 

4750  to  5160 

no  typical  value  -  user  must 

of  alternating  depth  and 
sound-speed  values  or  depth 
and  temperature  values 

m,  m/s 
ft,  °F 
m,  °C 

1440  to  1590 
30  to  85 

0  to  30 
at  the  surface 

enter  SSP  for  program  to  run 

BD,  ZBM 

Bottom  depth 

ft 

1  to 

largest  depth  in  SSP 

SD,  ZX 

Source  (transmitter)  depth 

ft 

1  to  1030 

20 

RD,  ZR 

Receiver  depth 

ft 

1  to  1030 

20 

TD,  ZTG 

Target  depth 

ft 

1  to  1030 

50 

FREQ,  F 

Frequency  of  transmitter; 

Hz 

75  to  105 

— 

center  frequency  of  receiver 
bandwidth 

kHz 

10'4  to  24.99 

3.5 

WAVE,  LWA 

Wave  height  (crest  to  trough) 

ft 

0  to  20 

4 

WIND,  VW1 

Wind  speed 

knots 

0.1  to  35 

15 

SLPB,  HS 

Signal  loss  per  bounce  at  the 
surface 

dB 

-2  to  20 

3,  calculated  internally  if 
negative 

LAT 

Latitude 

deg 

-90  to  90 

0 

BRLP,  AHB 

Bottom-reflection-loss  pro¬ 
file.  Array  of  alternating 
grazing-angle  and  bottom- 
reflection-loss  values 

deg,  dB 

not  checked 

zero  loss  at  any  grazing  angle 

XVRP,  BEAMX 

Transmit  vertical  response 
pattern  (i.e.,  beam  pattern) 

-  array  of  alternating  angle 
and  response  values,  where 
response  is  given  as  positive 
dB  down  from  the  maximum 
response,  defined  as  0  dB  at 
an  angle  of  0  deg 

deg,  dB 

not  checked 

omnidirectional  (i.e.,  0  dB 
for  all  angles);  if  the  array 
contains  one  element  only, 
the  response  pattern  sin(x)/x 
for  a  continuous  line  is  calcu¬ 
lated  (XVBM  needed) 

RVRP,  BEAMR 

Receive  vertical  response 
pattern 

same  as  for  XVRP 

(RVBM  needed  for  sin(x)/x 
response  pattern) 

XVBM,  DELPX 

Effective  vertical  beamwidth 
-  transmit 

deg 

0.1  to  180 

180 

RVBM,  DELPR 

Effective  vertical  beamwidth 
—  receive 

deg 

0.1  to  180 

180 

XDE,  PHDX 

Transmit  depression- 
elevation  angle  of  the 
main  axis  of  the  beam 
pattern 

deg 

-85  to  85 

5 

RDE,  PHDR 

Receive  depression  -  eleva¬ 
tion  angle  of  the  main  axis 
of  the  beam  pattern 

deg 

-85  to  85 

5 
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Table  B-l.  Continued 


ENVIRONMENTAL  PARAMETERS 

Name 

Description 

Units 

Range  Typical  Value;  Comments 

R,  RANGE 

Array  of  target  ranges  — 

kyd 

.001  to  1000  1  10  1  , 

may  be  entered  as  single 

values  or  in  groups  of  3: 

3rd  successive  value  may  be 
used  to  step  between  two 
preceding  values  - 

PLOT  Array  of  plot  flags  —  enter  —  —  all  flags  0 

as  PLOT  =  A  B  C 
A:  plot  propagation  loss 
=  1 :  incoherent  sum  and 
largest  arrival 
=  2 :  incoherent  sum  only 
=  3 :  largest  arrival  only 
B:  plot  reverberation 

=  1 :  pulse-averaged  reverb 
=  2 :  instantaneous  reverb 
=  3:  both  on  separate 
graphs 

C:  plot  signal  excess 
=  1  incoherent  sum  and 
largest  arrival 
=  2  incoherent  sum 
=  3  largest  arrival 

OUTP  =  1 :  limit  calculations  to  —  0  to  1  0 

propagation  loss  —  skip  re¬ 
verberation  and  detection 


REVERBERATION  PARAMETERS 

BBS,  MUB  Bottom-backscattering  dB/yd^//l  yd  -40  to -5  -27 

strength  if  internal  function 
is  used 

'y 

ABBS  Array  of  alternating  grazing  deg,  dB/yd  //I  yd  not  checked  array  all  zeroes 

angle  and  bottom-back- 
scattering  strength  values  — 
internal  function  not  used 

ASBS  Array  of  alternating  grazing  deg,  dB/yd^//l  yd  not  checked  array  all  zeroes 

angle  and  surface-back- 
scattering  strength  values  — 
do  not  use  if  internal  surface- 
backscattering  model  is 
desired 
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Table  B-l.  Continued 


REVERBERATION  PARAMETERS 


Name 

Description 

Units 

Range 

Typical  value;  Comments 

VBS,  MUV 

Volume-backscattering 
strength  of  the  water  column 
1  yd“  in  area  at  surface  ex¬ 
tending  to  depth  ZC 

dB/yd2//l  yd 

-300  to  0 

-49 

SLD,  ZC 

Scattering  layer  depth; 
maximum  depth  of  the 
water  column 

ft 

Oto  10,000 

2000 

SL 

Source  level 

dB//l  £(Pa//l  yd 

40  to  300 

240 

HBM,  DF.LTH 

Horizontal  beamwidth 

deg 

0.1  to  360 

360 

TBP 

Time  between  pings 

s 

1  to  200 

10 

PULS,  PULSE 

Pulse  length 

s 

104  to  100 

0.1 

BWR 

Bandwidth  of  receiver 

Hz 

104  to  1000 

100 

GDOP 

Doppler  gain 

dB 

-2  to  300 

-1,  internal  model  is  used  if 
GDOP  is  negative 

TCV,  TVC 

Target  closing  velocity  — 
used  only  with  internal 
doppler  model  when  GDOP 
=  -l 

Reverberation  output  flag 
=  1 :  print  special  reverber¬ 
ation  outputs 
=  2 ;  do  above  and  skip  de¬ 
tection  section 

knots 

-60  to  60 

0 

OUTR 

Oto  2 

0 

DETECTION  PARAMETERS 

DI 

Directivity  index 

dB 

Oto  70 

0 

NIN,  NSPEC 

Spectrum  level  noise  into 
the  receiving  array  (not 
needed  if  NOUT  is  used) 

dB/ 1  Hz//1  /U  Pa 

-80  to  100 

55 

NOUT,  NBWR 

Noise  in  the  receiver 
bandwidth  at  the  output 
of  the  beamformer  (super¬ 
cedes  NIN) 

dB// 1  /UPa 

-105  to  129 

65 

TGS 

Target  strength 

dB 

0  to  40 

15 

RDN 

Recognition  differential 

dB 

-20  to  20 

-2 

for  noise-limited  detection 
(corrected  by  program  for 
reverberation-limited 
detection) 
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Table  B-l.  Continued 


DETECTION  PARAMETERS 


Name 

Description 

Units 

Range 

Typical  Value;  Comments 

ARD 

Array  of  alternating  rever- 
beration-to-noise-ratio  and 
recognition-differential 
values 

dB,  dB 

not  checked 

RDN  =  -2 

ADT 

Array  of  alternating  rever- 
beration-to-noise-ratio  and 
detection-threshold  values 

dB,  dB 

not  checked 

RDN  =  -2 

CLIP 

=:1  reduces  signal  level 

2  dB  for  clipping  in  the 
signal  processor  (=  0,  no 
clipping  loss) 

0  to  1 

0 

SIG,  SIGMA 

Standard  deviation  of 
signal  differential  (only 
used  with  SP  =  5) 

dB 

0.5  to  6 

6 

PFA 

Probability  of  false 
alarm 

- 

10'12  to  0.1 

0.001 

SP 

Flag  to  denote  signal 
processing  model 

- 

1  to  5 

5 

=  1 :  linear  correlator,  known 
signal 

=  2 :  linear  correlator,  con¬ 
stant  signal 

=  3 :  quadrature  correlator, 
constant-amplitude  signal 
=  4:  quadrature  correlator, 
random-amplitude-and-phase 
signal 

=  5 :  energy  detector  for 
gaussian  signal 
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Table  B-2.  Data  Literals. 


♦ 

Table  B-2.  Data  Literals. 

Data  literal 

Resulting  control 

RUN 

Perform  calculations  for  the  values  of  the  current  input  parameters 
and  return  control  to  input  section. 

END 

End  program  execution. 

HEADER 

Use  the  next  card  image  for  output  headings. 

ABUT 

Abut  the  characters  in  the  next  card  image  to  the  current  charac¬ 
ters  in  the  heading. 

TTY 

Reduce  output  for  demand  terminal. 

NOTT* 

Print  normal  outputs,  defeat  TTY. 

NOPR 

Eliminate  most  printed  outputs. 

PRIN* 

Allow  printed  outputs,  defeat  NOPR. 

DR 

Print  detection  ranges. 

DRONLY 

Print  detection  ranges  only,  eliminate  other  outputs. 

NODR* 

Defeat  DR  and  DRONLY. 

GOPL 

Plot  previously  calculated  arrays  (appropriate  plot  flags  must  be  set). 

NOPL 

Allow  no  plotting  at  the  end  of  this  run,  clear  plot  flags. 

CW 

Label  the  plots  for  a  CW  signal. 

FM 

Label  the  plots  for  an  FM  signal. 

PRN* 

Label  the  plots  for  a  PRN  signal. 

LA* 

Use  the  propagation  loss  for  the  largest  arrival  to  compute 
detection  ranges. 

SUM 

Use  the  propagation  loss  for  the  incoherent  sum  to  compute 
detection  ranges. 

SAVY 

Write  outputs  in  Unit  1 1  (disk). 

SAVN* 

Do  not  write  outputs  on  a  disk. 

CAUN 

Use  standard  ray  tracing  without  caustic  corrections. 

CAUY* 

Use  caustic  corrections. 

XSUR  • 

Exclude  up-going  rays  at  the  source. 

XBOT 

Exclude  down-going  rays  at  the  source. 

RSUR 

Exclude  up-going  rays  at  the  receiver. 

RBOT 

Exclude  down-going  rays  at  the  receiver. 

TSUR 

Exclude  up-going  rays  at  the  target. 

TBOT 

Exclude  down-going  rays  at  the  target. 

The  program  is  initialized  to  this  condition  at  the  beginning  of  execution. 
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is  less  than  half  the  difference  of  the  first  two.  If  so,  the  third  element  is  used  as  a  stepping 
increment  between  the  first  two  elements.  Thus,  R  =  4  1  1  2  will  generate  the  sequence  4, 
6,  8,  10,  11.  Values  in  R  which  do  not  satisfy  the  stepping  principle  are  used  singly. 

Thus, 

R  =  1  5  1  4.5  7  10  2 
yields  the  sequence 

1,  2,  3,4,  4.5,  5,  7,  9,  10 
The  same  range  sequence  is  generated  by 
R  =  7  10  2  4.5  1  5  1. 

“7  10  2”  is  taken  as  a  stepping  sequence,  and  none  of  its  values  can  participate  in  another 
sequence,  such  as  “2  4.5  1 .”  A  maximum  of  90  values  can  be  used  in  R  and  the  maximum 
number  of  values  that  can  be  generated  for  range  in  LIRA  is  500.  For  plotting  purposes  it 
is  necessary  to  generate  range  values  closely  and  evenly  spaced. 

PLOT.  Array  of  plot  flags.  If  “PLOT”  is  entered  with  no  values  following,  then 
“PLOT  =  1  1  1”  is  assumed. 

PLOT  (1)  =  1 :  plot  propagation  losses  for  incoherent  sum  and  largest  arrival 

=  2:  plot  propagation  loss  for  incoherent  sum 
=  3  :  plot  propagation  loss  for  largest  arrival 
PLOT  (2)  =  1 :  plot  pulse-averaged  reverberatioif 

=  2:  plot  instantaneous  reverberation* 

=  3:  plot  both  pulse-averaged  and  instantaneous  reverberation  on 

separate  graphs 

PLOT  (3)  =  1 :  plot  signal  excess  for  the  incoherent  sum  and  the  largest  arrival 

on  the  same  graph 

=  2:  plot  signal  excess  for  the  incoherent  sum 
=  3:  plot  signal  excess  for  the  largest  arrival 

OUTR.  Reverberation  output  flag. 

=  0.  no  effect 

=  1 :  prints  bottom, surface  and  volume  reverberation  components,  total  contri¬ 
bution  from  each  ping  (if  more  than  one  ping)  and  total  reverberation  vs  target  range 
(and  corresponding  two-way  travel  time) 

=  2:  same  as  OUTR  =  1  but  detection  calculations  are  skipped. 

APT.  Detection  threshold  (DT)  vs  reverberation-to-noise  ratio  (RNR);  array  of 

Included  in  the  reverberation  plots  are  the  total  reverberation,  the  surface,  volume,  and  bottom 
components  of  reverberation,  and  the  ambient  noise  level. 


alternating  RNR  and  DT  values,  both  in  dB.  A  single  value  in  the  array  signifies  that  DT 
is  constant  at  that  value.  The  values  for  DT  are  measured  at  the  output  of  the  signal 
processor.  Detection  threshold  in  dB  is  the  ratio  of  the  signal-plus-noise  power  to  the 
noise  power  in  the  receiver  bandwidth  needed  to  achieve  a  specified  level  of  correct  de¬ 
cision  with  a  specified  probability  of  false  alarm  (PFA).  For  the  probability-of-detection 
models  in  LIRA  to  apply,  the  specified  level  of  correct  decision  is  probability  of  detection 
=  0.5  for  signal  excess  =  0  for  user-specified  PFA.  The  user  may  use  a  DT  corresponding  to 
another  probability  of  detection  (at  which  plotted  signal  excess  =  0),  but  the  printed 
values  for  probability  of  detection  vs  range  will  not  apply.  The  entering  of  the  ADT 
array  nullifies  previous  entries  of  either  RDN  or  ARD. 

ARP.  Recognition  differential  (RD)  vs  reverberation-to-noise  ratio  (RNR);  array 
of  alternating  RNR  and  RD  values,  both  in  dB.  A  single  value  in  the  array  signifies  that 
RD  is  constant  at  that  value.  The  values  for  RD  are  measured  at  the  output  of  the  beam- 
former.  RD  is  the  signal-to-masking-background  ratio  in  the  receiver  bandwidth  needed  to 
achieve  50%  probability  of  detection  for  a  given  probability  of  false  alarm.  ARD  is  a 
curve  representing  the  uniting  of  the  signal  processing  and  the  decision  processes  into  one. 

In  general  the  values  of  PFA  used  for  ARD  and  ADT  inputs  are  different  for  a  given  receiver 
because  the  signal  and  noise  distributions  are  different  before  and  after  the  signal  processor. 
Signal  excess  is  zero  for  50%  probability  of  detection.  Reverberation  plots  are  for  rever¬ 
beration  in  the  receiver  bandwidth  at  the  output  of  the  beam  former.  The  entering  of 
the  ARD  array  nullifies  previous  entry  of  either  RDN  or  ADT. 

Additional  Inputs 

Other  inputs  exert  specific  control  not  needed  by  the  general  user. 

NBBB,  REFLS.  Maximum  number  of  bottom  bounce  ray  cycles  for  which  ray 
tracing  calculations  will  be  performed.  If  this  parameter  is  not  entered,  a  maximum  of 
100  cycles  is  assumed,  but  calculations  are  limited  by  the  end  points  of  the  input  ranges 
to  increase  program  efficiency. 

ZONE,  ZONES.  Maximum  number  of  convergence  zone  ray  cycles  for  which 
ray  tracing  calculations  will  be  performed.  If  this  parameter  is  not  entered,  a  maximum 
of  100  cycles  is  assumed,  but  calculations  are  limited  by  the  end  points  of  the  input 
ranges  to  increase  program  efficiency. 

DUMP.  Flag  for  intermediate  printout. 

=  0:  no  effect 

=  1 :  lists  sound-speed  maxima  and  subsurface  ducts;  lists  sound  speeds 
used  for  ray  tracing  to  the  target;  lists  surface  duct  propagation  loss,  surface  and  volume 
reverberation,  and  two-way  travel  time  to  the  target. 

FLAG.  Array  containing  flags  which  control  diagnostic  printouts. 

FLAG  (N)  =0,  N=  1,..,  6:  no  effect. 

FLAG  (1)  applies  to  caustics  generated  in  both  the  propagation  loss  and  rever¬ 
beration  calculations. 

=  1:  print  parameters  for  each  caustic 
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=  2  :  add  calculations  for  each  iteration  for  each  caustic 

=  3  :  add  details  for  each  iteration 

FLAG  (2)  applies  to  ray  interpolations  for  propagation  loss  to  the  target. 

=  1  :  print  parameters  for  individual  rays  for  both  bottom  bounce 
and  convergence  zone  paths 

=  2 :  add  the  interpolations  for  range  after  each  successive  ray 

=  -1 :  only  print  backscattering  strengths  for  each  ray 

=  -2:  print  a  dump  of  array  “ARP,”  the  array  containing  the  ray 
segment  parameters  for  all  rays  (all  bottom  bounce  and  convergence  zone  ray  segments 
for  bottom  surface  and  volume  backscatterers  and  for  the  target) 

FLAG  (3)  >\  :  allow  FLAG  (1 )  to  apply  to  reverberation  caustics. 

>2  :  add  printing  of  minimum  and  maximum  two-way  travel  times 
for  each  family  of  rays,  ray  cycle,  and  ray  type  for  each  backscatterer. 

>2 :  add  printing  of  lengthy  reverberation  interpolations  as  they  are 
being  calculated  vs  time. 

FLAG  (4)  =  1  :  print  a  dump  of  arrays  “ABB”  and  “ACZ,”  the  arrays  containing 

the  bottom  bounce  and  convergence  zone  vertex  sound  speeds  used  in  ray  tracing. 

FLAG  (5)  limits  calculations  to  a  particular  backscatterer;  all  other  calculations 
are  shunted. 

=  1  :  bottom 

=  2 :  bottom  bounce  paths  to  the  surface 

=  3 :  refractive  paths  to  the  surface 

=  4 :  bottom  bounce  paths  to  the  volume  or  target 

=  5  :  refractive  paths  to  the  volume  or  target 

FLAG  (6)  =  K ;  K  =  1 ,  2,  3,  4.  Limit  calculations  to  ray  type  “K”  for  both 

propagation  loss  to  the  target  and  reverberation. 
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APPENDIX  C 

INPUT  FORMAT 

The  LIRA  program  accepts  all  inputs  through  SUBROUTINE  READIN.  READIN 
contains  only  two  read  statements.  The  input  characters  are  scanned  in  a  manner  similar 
to  NAMELIST  to  determine  if  a  string  of  characters  is  an  input  name  or  an  input  para¬ 
meter  value  (a  number).  READIN  does  not  use  NAMELIST  for  the  following  reasons  : 

•  It  is  important  to  keep  track  of  the  variables  which  are  changed  from  run  to 
run  in  order  to  reduce  redundant  calculations. 

•  Column  1  on  the  input  record  can  now  be  used. 

•  Data  literals  which  have  no  numerical  value  can  be  entered. 

The  user  enters  single-value  variables  by  typing  the  name  of  the  variable,  followed 
by  a  space  or  an  equal  sign  “=”  and  then  by  the  numerical  value  in  F,  I,  E,  or  D  format. 

All  input  values  are  converted  to  single-precision  floating  point  numbers.  There  can  be  no 
imbedded  blanks  in  either  the  name  or  value,  although  there  may  be  more  than  one  blank 
separating  the  name  and  value.  A  name  and  its  value  must  be  on  the  same  line.  If  more 
than  one  name-value  set  occurs  on  a  line,  the  sets  are  to  be  separated  by  a  comma  or 
a  slash  “/  .  There  is  no  need  to  use  a  comma  or  a  slash  as  a  separator  between  sets  occur¬ 
ring  on  different  lines.  The  parameters  can  be  input  in  any  order  or  not  input  at  all,  as  the 
user  desires. 

The  user  enters  arrays  in  a  fashion  similar  to  that  used  for  single-value  variables, 
by  first  typing  the  name,  then  a  space  or  equal  sign  “=”,  and  then  the  array  values  in  F, 

I,  E,  or  D  format.  Successive  array  values  are  separated  by  one  or  more  spaces.  If  the 
array  is  too  long  for  one  line,  as  many  lines  as  are  required  can  be  used  by  typing  the 
continuation  symbol  “*  ”,  that  is,  an  asterisk  and  a  space,  as  the  first  characters  on  the 
successive  lines.  An  array  is  terminated  by  a  comma  or  a  slash  or  by  starting  a  new  line 
without  the  continuation  symbol. 

Data  literals  are  names  which  are  not  followed  by  values.  To  enter  the  name  is  to 
set  a  flag  or  to  follow  some  instructions  before  reading  the  next  input. 

Figure  C-l  illustrates  some  of  the  possible  ways  data  may  be  entered.  Column 
1  is  indicated  by  the  arrow.  The  line  numbers  are  merely  for  reference  in  this  discussion. 

Line  1  in  Fig.  C-l  shows  several  examples  of  separators  and  spacings  and  the  use  and  non¬ 
use  of  the  equal  sign.  Line  2  contains  the  data  literal  “DR,”  which  will  cause  a  flag  to  be 
set  so  that  detection  ranges  will  be  calculated.  Line  3  is  another  data  literal,  “HEADER,” 
which  causes  the  next  card  to  be  read  as  heading  information.  Line  5  shows  how  an  array 
of  target  ranges  might  be  entered.  The  user  need  not  specify  the  number  of  values  in  an 
array;  the  program  automatically  counts  them.  Lines  7-10  illustrate  the  entering  of  values 
in  array  “AHB"’,  employing  continuation  cards.  The  latest  entry  of  any  parameter  supersedes 
any  value  it  may  have  been  assigned  by  previous  entries.  On  card  1 1  “TGS  =  1 5”  is  super¬ 
seded  by  “TGS  =  20.”  On  Line  1 2  PFA  is  assigned  value  in  E  format  and  on  Line  13  in 
D  format.  The  data  literal  on  Line  12  (“RUN”)  will  cause  immediate  exit  from  the  READIN 
subroutine,  calculations  will  be  performed  with  the  current  data  set,  and  control  will  be 
returned  to  READIN  to  read  the  new  value  for  PFA  on  Line  13.  The  “RUN”  on  Line  14 
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will  again  cause  calculations  to  be  performed,  this  time  only  in  the  detection  section,  since 
PFA  is  a  detection  parameter.  Line  1 5  ends  the  program  with  “END,”  and  control  is 
returned  to  the  executive. 

Col.  1 

i 

1  FREQ=1„  ZX=20/ZTG  300.,PHDX=5  ,PHDR  .5 

2  DR 

3  HEADER 

4  THIS  IS  A  HEADING 

5  R  1  10  1  10  20  2  19  17  11  15  2 

6  SSP  0  5000  100  5001.8  3000  4950  10000  5050 

7  AHB  0  5  5.0  8.5  10.  12 

8  *  20  13  40  14 

9  *55  15 

10  *90  15 

11  TGS=15,  TTY,  TGS=20 

12  PFA=1.E-5,  RUN 

13  PFA=1  .D-4 

14  RUN 

15  END 


Figure  C-l.  Example  of  some  options  in  the  input 
format. 


A  data  set  equivalent  to  that  represented  in  Fig.  C-l  is  shown  in  Fig.  C-2. 


Col.  1 
i 

HEADER 

THIS  IS  A  HEADING 
FREQ=1,  ZX=20,  ZTG=300,  TGS=15 
PFA=. 00001,  PHDX=5,  PHDR=-5 
DR,  TTY 
R  =  1  20  1 

SSP=0  5000  100  5001.8  3000  4950  10000  5050 
AHB=0  5  5  8.5  10  12  20  13  40  14  55  15  90  15 
RUN 

PFA=.0001,  RUN 
END 

Figure  C-2.  Rewritten  form  of  the  data  in  Figure  C-l. 
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APPENDIX  D 


r 


COMPARING  LIRA  WITH  OTHER  MODELS 

The  most  sensitive  aspect  of  active-sonar  performance  prediction  is  the  calculation 
of  propagation  loss.  Predictions  are  equally  sensitive  to  underlying  theory  and  algorithms 
implementing  the  theory.  Within  any  practical  application  (i.e.,  computer  program)  of  a 
theoretical  approach  for  calculating  propagation  loss,  a  host  of  idiosyncratic  cases  exist. 
Each  case  requires  individual  treatment  and  numerical  approximation  when  no  analytical 
solution  exists  (usually  the  case).  Because  of  the  great  complexity  of  the  calculations,  some 
of  the  idiosyncracies  inherent  in  the  method  are  aften  overlooked.  Comparison  of  a  pro¬ 
pagation  loss  program's  outputs  with  data  often  looks  good  when  the  variance  of  the  data 
is  considered.  A  program  is  sometimes  rejected  when  it  fails  to  agree  with  data  and  vice 
versa 

There  is  no  simple  way  of  validating  a  program.  Perhaps  the  most  significant 
factors  in  the  validating  of  a  propagation  loss  program  are 

•  Recognition  of  the  validity  of  the  theoretical  approach. 

•  Lack  of  obvious  spurious  predictions. 

•  Good  comparison  with  data. 

•  Wide  range  of  use  in  the  acoustic  propagation  community. 

•  Agreement  with  other  validated  programs. 

The  FACT*  program  easily  falls  into  the  category  of  a  validated  program.  It  has 
its  shortcomings,  as  does  every  propagation  loss  program,  but  it  is  a  useful  and  much  used 
tool  for  passive-sonar  applications. 

FACT  (Fast  Asymptotic  Coherent  Transmission  model),  developed  by  the  former  AESD  (Acoustic 
Environmental  Support  Detachment)  at  Office  of  Naval  Research  (Ref.  9),  computes  incoherent  and 
semi-coherent  propagation  losses  using  ray  theory  with  corrections  for  both  smooth  and  cusped 
caustics.  FACT  is  the  “Navy  Standard”  model  for  lay  models.  FACT  represents  the  sound-speed 
profile  by  straight  lines  connecting  the  points.  Because  of  the  discontinuous  gradients  at  the  profile- 
segment  interfaces,  numerous  false  caustics  would  be  generated  by  standard  ray  theory.  Instead,  for 
a  particular  ray  family,  FACT  fits  range  vs  starting  angle  (R-0points)  with  a  curve  (R-0curve)that  is 
limited  to,  at  most,  one  extremum.  The  relatively  small  extrema  in  the  R-0  curves  caused  by  the 
gradient  discontinuities  are  washed  out  by  the  dominating  effect  of  the  real  caustic.  FACT  therefore 
forces  the  condition  of,  at  most,  one  caustic  per  ray  family.  Other  real  caustics  in  a  ray  family  are 
truly  obliterated  because  FACT  uses  the  fitted  R-0  curve  to  calculate  propagation  loss  rather  than 
the  R-0  points.  FACT  is  perhaps  the  only  ray  program  to  calculate  the  acoustic  field  at  and  near 
horizontal  cusps.  FACT  uses  a  surface  duct  model  developed  for  FNWC  (Ref.  10).  For  receivers 
below  the  surface  layer,  FACT  uses  the  surface  duct  propagation  loss  plus  a  constant.  When  source 
and  receiver  are  located  near  the  deep  sound  channel  axis,  FACT  uses  a  wave  correction  for  the 
rays  which  obtain.  The  maximum  angle  for  ray  tracing  in  FACT  is  30  deg:  thereafter  an  approximation 
method  is  used  for  propagation  loss.  Only  the  internal  bottom  loss  functions  may  be  used  in  FACT, 
there  is  no  provision  for  entering  a  new  one.  Neither  can  a  beam  pattern  be  entered. 
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Figure  D-l  (a).  Sound-speed  profiles  used  to  compare  LIRA,  PLRAY,  FACT  and  AP2  (from  Ref.  1 1). 
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Recently,  Bartberger  at  NADC  compared  his  newly  developed  ray  theory  program, 
PLRAY,  with  FACT  and  a  normal  mode  program,  AP2.t  He  compared  the  semi- 
coherent  propagation  losses  from  PLRAY  and  FACT  with  each  other  and  with  smoothed 
coherent  propagation  loss  from  AP2.  He  also  compared  incoherent  propagation  loss 
curves  from  PLRAY  and  FACT. 

The  method  used  here  to  partially  validate  LIRA  is  to  compare  LIRA’S  incoherent 
propagation  loss  predictions  with  those  of  PLRAY  and  FACT  for  the  14  cases  Bartberger 
used.  Bartberger’s  normal  mode  results  are  used  in  several  cases  to  illustrate  certain  salient 
aspects  of  the  comparison.  The  sound-speed  profiles  used  by  Bartberger  are  re-plotted  in 
Fig.  D-l.  Whenever  “NADC  SSP  N”  is  referred  to  from  now  on,  the  reader  will  find  the 
appropriate  plot  in  Fig.  D-l,  where  “N”  is  the  number  of  the  sound-speed  profile.  The 
bottom  loss  functions  of  Ref.  10  are  not  reproduced  here;  bottom  loss  values  in  the 
LIRA  runs  are  found  in  Table  D-l.  Figures  D-2  to  D-l 6  are  duplications  of  Bartberger’s 
propagation  loss  plots,  with  plots  of  LIRA’S  propagation  loss  curves  superposed  appro¬ 
priately  for  comparison.  The  propagation  losses  from  LIRA  are  interpolated  for  ranees 
from  0.2  kyd  to  100  kyd  in  steps  of  0.2  kyd.  It  is  not  known  what  interpolation  interval 
was  used  for  PLRAY  and  FACT. 

The  14  cases  will  now  be  examined  one  at  a  time,  with  comments  about  the 
different  models.  It  is  hoped  that  the  theories  or  algorithms  of  these  models  have  not 
been  misinterpreted  in  the  course  of  attempting  to  explain  the  differences  between  them 
and  LIRA.  The  bias  of  this  report  is  to  validate  LIRA,  and  the  comments  herein  regarding 
the  14  cases  are  not  intended  as  a  repudiation  of  either  of  the  other  fine  ray  models. 

CASE  1. 

NADC  SSP  1  is  used  with  source  depth  at  80  ft,  target  depth  at  100  ft,  and  frequency 
at  105  Hz.  Incoherent  propagation  loss  is  presented  for  LIRA,  PLRAY,  and  FACT  in 
Fig.  D-2. 


PLRAY  (Ref.  1 1 ),  developed  by  Charles  Bartberger  at  NADC  (Naval  Air  Systems  Command),  calculates 
incoherent  and  semi-coherent  propagation  losses  using  ray  tracing  with  corrections  for  smooth  caustics. 

The  sound-speed  profile  is  fitted  with  curvilinear  segments  having  continuous  gradients  at  the  interfaces. 
PLRAY  searches  for  caustics  in  ray  sub-families  defined  both  by  limiting  rays  and  certain  gradient 
changes.  PLRAY  finds  and  interpolates  a  maximum  of  one  caustic  per  ray  sub-family  even  though 
more  caustics  may  exist.  The  caustic  chosen  by  PLRAY  has  the  dominant  field.  The  range  regions 
of  the  other  caustics  are  interpolated  by  standard  methods.  PLRAY  merges  caustic  interpolation 
smoothly  with  ray  interpolation.  Horizontal  cusped  caustics  are  treated  by  offsetting  source  and 
receiver  in  depth  slightly  so  that  the  cusps  cannot  occur.  Normal  mode  methods  are  used  for  propa¬ 
gation  loss  in  the  surface  duct.  PLRAY  has  a  leaky  duct  model  for  transition  between  good  and  no 
propagation  in  the  surface  duct.  PLRAY  has  an  internal  set  of  bottom  loss  curves,  which  it  inter¬ 
polates  for  frequency.  Bottom  loss  may  also  be  read  in.  PLRAY  accepts  beam-pattern  inputs. 

tAP2  is  a  normal  mode  program  newly  developed  at  NADC,  still  undocumented  (Ref.  II).  AP2  uses 
a  maximum  of  500  modes  to  calculate  coherent  propagation  loss.  For  comparison  with  semi-coherent 
models,  AP2’s  curve  of  intensity  vs  range  is  smoothed  with  a  weighted  sliding  window  spanning  nine 
terms.  AP2  is  compatible  with  ray  programs  in  that  it  can  take  as  inputs  arbitrary  sound-speed  profiles 
and  bottom  loss  functions.  AP2  was  compared  with  other  normal  mode  programs  and  gave  identical  results 
for  the  same  cases. 
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HE aoer 
NADC  SSP  1 

SSP  0  5030.55  233.5  5033.63  395.6  5221.67  433.4  5021.44  597.1  5016.64 

*  755.9  5207.44  777. v  5^05.64  1096.4  50G1.75  1214.8  4999.05 

*  1520.7  4992.35  3005.3  4V2G.79  4011.7  4693.38  5014.1  4698.32 

*  8027.9  4934.30  10200  496C  ItOCO  4y90  14000  5020  15000  5035 
HS  C 

OOTP  1 
R  .2  100  .2 
FREG  105 

AHB  0  0  5  C  15  2  25  2  35  15  45  17  65  13  75  12  90  11.5 
2  TO  80,2*  100, RoN 
HEADER 
NADC  SSP  2 

SSP  0  5045  2CC  5025  400  4975  oOO  49o5  SO 0  4960  UOO  4943  2000  4932 

*  2400  4916  3CCC  4902  3400  4897  4000  4896  46Cu  4900  6CC0  4914 

*  8000  494 C  14000  5030  1800C  5090 
FREG  100 

AHB  0  2  10  2  2C  6  30  7.5  40  10  55  12  50  12 

7TG  300, ZX  60, RoN 

ZX  200, RUN 

HEADER 

NADC  SSP  3 

SSP  0  5005.9  32.8  5CG4.91  65.6  5005.9  9a. 4  4990.15  164.0  4948.81 

*  246.1  4933.39  328.1  4921.91  410.1  4914.69  492.1  4909.44 

*  964.2  4652.05  1640.4  4359.57  1968.5  4852.35  2256.6  4850. G6 

*  2624.7  4650.71  360t.9  4857.93  4921.2  4872.04  5741.5  4881.55 

*  6561.7  4692.05  6202.1  4916.33  9842.5  4942.56  13123.3  4999.33 

*  16404.2  5055.04  17762.4  5084.96 

AHB  0  0  10  0  20  3  35  7  45  3  50  10  90  10 
7  TG  300, ZX  60, RUN 
ZX  1000, PUN 
FREG  300 

AHB  0  2  10  2  12  2.5  20  4.5  35  7.5  45  9.5  55  11  90  11 
RUN 

HEADER 
NADC  SSP  4 

SSP  0  5019  275  5025  s50  5G0u  1300  4680  2000  4fc55  4000  4865  6000  4885 

*  1C000  4943  16000  5G50 
FREG  ISO 

AHB  0  0  20  0  25  2  30  4  55  8  9o  6 

7TG  200, ZX  60, PON 

ZX  300, RUN 

ZX  1000, RUN 

HEADER 

NADC  SSP  5 

SSP  0  4921.91  32.8  4915.67  65. t  4911.0?  98.4  4911.0a  160.8  4910.1 

*  242.8  4906.13  321.5  4694. u2  465.6  4e92.3e  6<.6.3  4o93.e9  967.8  4894.68 

*  1607.6  4899.6  1932.4  4697.63  3231.6  4*84.18  3681.2  4680.57 

*  4652.4  4668.1  1  6456.7  491  l.  .  75  8067.6  4933.06  9005.9  4944.54 
FREG. 100 

AHB  0  C  10  0  20  3  35  7  45  8  50  1C  90  10 
Z  TG  300. ZX  100, RUN 
Zx  1000, RON 


Table  D-l .  LIRA’S  inputs  for  the  14  cases  in  this  Appendix. 
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PEACE* 

NACC  SSP  6 

SEP  C  491C.8  98.4  4V1U.7  12-.7  4869.75  377.3  4827.6  964.2  4839.9 

*  14  32.7  4fc  32. 3l  492  1.3  4667.2  9842.3  4443.6  17782.2  50o3.7 
f  8  E  Q  20 

A  HD  C  6  90  6 

2TG  300, 7x  ICCO.RUN 

HEA0E8 

NACC  SSP  7 

SSP  0  5062. yO  32.6  5063.62  yO  5022.76  96.4  5023.62  246.1  4079.33 
«  360.9  4977.69  574.2  4976.67  1146.3  4976.05  1640.4  4979 

*  2952.6  4995.41  5000  5029.18  7CCC  5062.65  9842.5  5111.88 
f«E«.  100 

AND  0  0  10  0  20  3  35  7  45  6  5,  1C  90  1C 
7TG  3 JC i 7  X  100, HUN 
ZX  100C.RUN 
END 


Table  D-l.  Continued. 
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The  three  programs  agree  in  the  bottom  bounce  region  (4<  R  <  75  kyd  and  R>85 
kyd).  There  are,  however,  subtle  differences  in  the  direct-path  and  convergence  zone 
regions. 

All  three  models  use  their  own  specialized  surface  duct  models  derived  from  data 
and  normal  mode  theory.  The  surface  duct  propagation  losses  for  ranges  less  than  4  kyd 
are  slightly  different  for  each  model. 

At  75  kyd  there  appears  to  be  the  onset  of  a  convergence  zone.  Actually,  the 
decrease  in  propagation  loss  is  produced  by  rays  which  converge  but  do  not  cross;  that  is, 
there  are  no  caustics  at  all  produced  by  this  geometry,  at  least  for  LIRA.  The  shape  of 
the  FACT  curve  looks  like  part  of  an  Airy  function  curve  that  has  been  cut  off.  It  is 
likely  that  the  FACT  program  established  the  presence  of  a  reversed  caustic  (i.e.,  range 
vs  starting  angle  has  a  maximum)  and  interpolated  for  that  caustic  on  the  illuminated 
side  for  only  those  ranges  which  had  rays  arriving.  The  shadow  zone  extends  past  50  kyd, 
and  again  rays  contributing  to  the  field  are  excluded  because  they  fall  within  the  boundaries 
of  the  caustic  field. 

In  contrast,  LIRA  and  PLRAY  have  not  established  the  presence  of  a  caustic  and 
use  standard  ray  tracing  for  the  acoustic  field.  For  this  reason  the  tail  of  the  zone  extends 
farther  than  FACT’S. 

It  is  not  the  purpose  of  this  discussion  to  make  a  choice  between  models.  Because 
of  the  different  representations  of  the  sound-speed  profile,  one  model  can  certainly 
establish  the  presence  of  a  caustic  and  others  not.  Indeed,  there  is  cause  for  encouragement 
that  the  fundamentally  different  algorithms  agree  so  well  in  their  predictions. 

CASE  2. 

NADC  SSP  2  is  used  with  source  depth  at  60  ft,  target  depth  at  300  ft  and  fre¬ 
quency  at  100  Hz.  Incoherent  propagation  loss  for  LIRA,  PLRAY,  and  FACT  is  presented 
in  Fig.  D-3. 

There  is  good  agreement  between  the  three  models  in  the  direct  path  and  bottom 
bounce  regions  as  well  as  the  leading  edge  (shadow  side)  of  the  convergence  zone.  Slight 
variations  in  sound-speed  profile  representation  can  cause  the  shifting  of  the  intensity  of 
the  acoustic  pressure  field  in  the  caustic  regions.  Both  PLRAY  and  FACT  use  approxima 
tions  for  the  position  of  the  caustics  rather  than  the  iterative  method  of  LIRA,  probably 
causing  slight  changes  in  the  positions  and  shapes  of  the  several  caustic  fields  which  make 
up  the  composite  field.  The  average  magnitudes  are  the  same  over  the  region  of  the  con¬ 
vergence  zone  from  68  to  74  kyd. 

All  three  models  forbid  interpolation  of  a  caustic  on  the  illuminated  side  for  those 
ranges  not  reached  by  rays  in  the  same  ray  family  from  which  the  caustic  is  derived. 

This  results  in  the  sharp  increase  in  propagation  loss  at  73-74  kyd  for  all  three  programs. 

For  LIRA  and  PLRAY  two  caustics  were  truncated  in  this  manner,  resulting  in  the  small 
notch  at  74  kyd.  FACT  truncated  one  caustic  at  the  tail  of  the  convergence  zone.  Perhaps 
if  FACT  had  determined  the  other  caustic  at  the  tail  of  the  zone,  the  shape  of  the  zone 
would  have  been  more  like  PLRAY’s. 
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CASE  3. 


The  incoherent  propagation  losses  in  Fig.  D-4  are  for  the  same  sound-speed  profile 
as  CASE  2;  source  depth  is  changed  to  200  ft.  The  same  remarks  apply,  approxi¬ 
mately,  to  this  case  as  CASE  2,  since  there  are  no  novel  differences. 

CASE  4. 

NADC  SSP  3  is  used  with  source  depth  at  60  ft,  target  depth  at  300  ft  and  fre¬ 
quency  at  100  Hz.  Incoherent  propagation  loss  for  LIRA,  PLRAY,  and  FACT  is  presented 
in  Fig.  D-5. 

The  small  peak  at  4  kyd  in  LIRA  is  not  a  caustic,  just  a  narrow  bundle  of  rays. 

No  novel  differences  occur  at  convergence  zone  ranges. 

Neither  PLRAY  nor  FACT  traces  bottom  bounce  rays  steeper  than  30  deg.  They 
both  use  the  same  approximation  to  compute  propagation  loss  for  ranges  corresponding 
to  bottom  bounce  with  steeper  rays.  LIRA  traces  rays  at  angles  up  to  89  deg,  nearly 
vertical.  This  can  account  for  the  slight  discrepancy  in  bottom  bounce  propagation  loss 
between  LIRA  and  the  other  two  models  for  ranges  between  5  and  10  kyd. 

CASE  5. 

Incoherent  propagation  loss  for  NADC  SSP  3,  for  a  frequency  of  100  Hz,  source 
depth  at  1000  ft  and  target  depth  at  300  ft,  is  presented  in  Fig.  D-6  for  LIRA,  PLRAY, 
and  FACT.  Smoothed  coherent  propagation  loss  for  AP2  is  presented  in  Fig.  D-7. 

The  most  noteworthy  differences  between  LIRA,  PLRAY,  and  FACT  are  the 
heights,  widths,  and  numbers  of  the  caustic  peaks  in  the  convergence  zone  region.  PLRAY’s 
leftmost  peak  is  narrower,  higher,  and  slightly  to  the  left  of  the  corresponding  peak  in 
LIRA  and  FACT.  FACT’S  leftmost  peak  follows  the  opposite  trend,  leaving  LIRA’S 
leftmost  peak  an  intermediate  between  those  of  PLRAY  and  FACT.  A  comparison 
with  the  normal  mode  curve  for  CASE  5  in  Fig.  D-7a  shows  that  LIRA’S  prediction  con¬ 
forms  better  to  the  AP2’s  prediction  than  those  of  PLRAY  and  FACT.  The  same  is  true 
of  the  caustic  peak  at  59  kyd  in  the  LIRA  curve. 

The  peak  at  56  kyd  between  the  two  most  prominent  peaks  in  the  LIRA  curve 
does  not  have  any  counterpart  in  the  normal  mode  curve.  It  may  be  that  secondary  caustics 
in  a  ray  family  which  occur  because  of  slight  sound-speed  profile  nuances  should  be  excluded 
from  interpolation.  In  PLRAY  the  secondary  caustic  occurs  at  54  kyd.  From  its  shape  it 
appears  to  have  escaped  the  caustic  search.  The  small  peak  was  probably  obtained  by  ray 
interpolation.  FACT’S  treatment  of  this  secondary  caustic  is  consistent  with  the  theory  that 
more  than  one  caustic  exists  in  this  ray  family.  Assuming  that  the  secondary  caustic  does 
occur  in  the  ray  family,  FACT  would  fit  a  curve  with  a  single  maximum  through  the  R-0 
points  containing  both  the  primary  and  secondary  maxima.  FACT  then  would  use  the 
fitted  curve  to  calculate  the  parameters  of  a  single  caustic,  resulting  in  a  wider,  lower  peak 
than  if  the  secondary  caustic  did  not  exist. 

CASE  6. 

The  environment  and  geometry  are  the  same  as  the  previous  case.  The  only 
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difference  is  that  frequency  is  300  Hz.  The  incoherent  propagation  loss  curves  are  shown 
in  Fig.  D-8.  Bartberger  evidently  included  this  case  in  his  comparisons  to  show  that  FACT 
does  not  interpolate  its  internal  bottom-loss  curves  for  frequency. 

In  the  direct-path  region,  LIRA  resolves  two  caustics;  PLRAY  and  FACT  each 
resolve  one.  Perhaps  this  is  another  case  of  more  than  one  caustic  per  ray  family.  At 
300  Hz  the  caustics  in  the  convergence  zone  region  are  more  highly  resolved,  lending 
visual  support  to  the  theory  of  CASE  5.  The  flattened  regions  following  each  of  FACT’S 
two  main  caustics  in  the  convergence  zones  are  probably  not  obtained  by  ray  interpolation 
but  by  an  approximation  for  the  “beat”  region  of  the  Airy  function,  i.e.,  the  illuminated 
side  of  the  caustic  outside  of  the  boundary  layer.  (See  Fig.  4.)  Neither  LIRA  nor  PLRAY 
uses  this  approximation.  Rather,  they  interpolate  ray  contributions  immediately  outside 
the  caustic  boundary  layer.  It  is  assumed  that  FACT  interpolates  its  approximation  for 
the  “beat”  region  over  the  range  interval  for  which  rays  arrive  from  the  ray  family  con¬ 
taining  the  caustic. 

Bartberger  was  unable  to  present  a  normal  mode  comparison  for  this  frequency 
because  of  the  excessive  number  of  modes  required. 

It  is  not  understood  why  PLRAY’s  caustic  at  58  kyd  is  not  prominent. 

CASE  7. 

NADC  SSP  7  is  used  with  source  depth  at  60  ft,  target  depth  at  200  ft,  and  fre¬ 
quency  at  1 50  Hz.  The  incoherent  propagation  loss  curves  are  shown  in  Fig.  D-9. 

Both  the  source  and  the  target  are  in  the  surface  duct,  so  all  of  the  convergence 
zone  caustics  are  formed  by  rays  which  (eventually)  strike  the  surface.  For  this  zone 
there  is  a  possibility  of  four  caustics  occurring,  but  because  of  the  close  proximity  of  the 
source  and  target  to  the  surface,  they  are  not  all  resolved. 

As  in  previous  cases  LIRA  resolves  more  caustics  than  PLRAY  or  FACT.  There 
is  a  hint  of  a  second  caustic  in  the  PLRAY  curve.  In  the  normal  mode  run,  not  shown, 
two  caustics  are  resolved,  with  the  second  more  predominant  than  the  First. 

As  in  previous  cases,  the  onset  of  the  zone  (leading  edge  or  shadow  region)  is 
steeper  for  PLRAY  than  it  is  for  LIRA,  FACT,  and  AP2. 

CASE  8. 

The  only  difference  between  this  and  the  previous  case  is  that  target  depth  is  300 
ft,  located  below  the  duct.  The  propagation  loss  curves  are  shown  in  Fig.  D-10. 

LIRA  and  PLRAY  both  show  direct-path  propagation  at  longer  ranges  than 

FACT. 

LIRA  shows  three  caustics  fully  resolved  in  the  convergence  zone  and  the  hint  of 
a  fourth.  PLRAY  shows  three  caustics  and  FACT  only  two.  Again,  the  onset  of  PLRAY’s 
leading  caustic  is  steeper  than  that  of  either  LIRA  or  FACT.  The  positions  of  PLRAY's 
caustic  peaks  are  at  least  1  kyd  less  than  corresponding  peaks  in  LIRA  and  FACT. 
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CASE  9. 


This  case  is  the  same  as  the  previous  one,  but  with  the  target  depth  at  1000  ft. 

The  incoherent  propagation  loss  curves  are  shown  in  Fig.  D-l  1. 

Each  of  the  caustics  for  the  four  ray  types  is  fully  resolved  in  both  LIRA  and 
PLRAY.  It  is  not  understood  why  FACT  does  not  resolve  four  caustics,  since  each  caustic 
belongs  to  a  different  ray  type.  Perhaps  the  algorithm  in  FACT’S  semi-coherent  model 
used  to  combine  ray  types  when  source  or  receiver  is  near  the  surface  was  applied  here. 

The  peak  at  40  kyd  in  PLRAY’s  curve  appears  to  be  spurious. 

CASE  10. 

NADC  SSP  5  is  used  with  source  depth  at  100  ft,  target  depth  at  300  ft,  and  fre¬ 
quency  at  100  Hz.  Incoherent  propagation  loss  is  presented  in  Fig.  D-12. 

LIRA  and  PLRAY  agree  in  the  direct-path  region  of  propagation  loss.  FACT, 
however,  is  closer  to  the  normal  mode  prediction  (not  shown). 

The  three  models  differ  in  the  convergence  zone,  although  LIRA  and  PLRAY  are 
the  two  models  in  closest  correspondence.  This  convergence  zone  is  a  complex  super¬ 
position  of  eight  or  more  caustics  derived  from  totally  refracted  rays  and  surface-reflected 
rays.  Mere  observation  of  the  propagation  loss  curve  yields  no  new  information  about  the 
structure  of  the  model. 

CASE  11. 

This  is  the  same  as  the  previous  case,  with  source  depth  moved  to  1000  ft.  The 
propagation  loss  curves  are  presented  in  Fig.  D-l 3. 

In  this  example  propagation  occurs  in  a  submerged  duct  as  well  as  the  deep  paths. 
LIRA  and  PLRAY  both  Find  caustics  which  are  nearly  horizontal  in  the  submerged  duct 
and  truncate  their  Fields  on  the  illuminated  sides  at  ranges  for  which  there  are  no  ray 
arrivals.  Caustics  for  which  the  illuminated  region  is  both  increasing  and  decreasing  with 
range  are  found  by  LIRA  and  PLRAY.  As  a  result  the  propagation  loss  curves  from  both 
LIRA  and  PLRAY  exhibit  a  curious  scalloped  shape. 

There  are  several  possibilities  why  FACT’S  curve  does  not  show  the  same  pro¬ 
perties: 

•  FACT  may  have  some  limit  at  which  it  rejects  contributions  from  nearly 
horizontal  caustics. 

•  FACT  may  not  have  traced  rays  in  the  duct  because  of  insufficient  duct 
strength. 

•  FACT  may  have  used  its  wave  approximation  for  a  deep  duct. 

Whatever  the  approach  FACT  used,  its  curve  matches  the  normal  mode  curve  (not 
shown)  better  than  LIRA  or  PLRAY. 

The  peak  at  21  kyd  in  the  curve  for  LIRA  is  due  to  rays  that  converge  but  do 
not  cross.  No  caustic  solution  is  obtainable. 
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CASE  12. 


NADC  SSP  6  is  used  with  source  depth  at  1000  ft,  target  depth  at  300  ft,  and 
frequency  at  50  Hz.  Incoherent  propagation  loss  curves  for  LIRA,  PLRAY,  and  FACT  are 
presented  in  Fig.  D-14.  A  smoothed  coherent  propagation  loss  curve  for  AP2  is  presented 
in  Fig.  D-7b. 

LIRA  computes  in  the  direct-path  region  a  strong  caustic  that  does  not  occur  in  the 
other  three  programs.  PLRAY’s  direct  path  seems  to  cut  off  early  in  range. 

LIRA  and  PLRAY  compare  better  with  AP2  than  FACT  in  the  convergence  zone 
regions.  All  the  predominant  caustics  are  nearly  horizontal.  It  appears  that  FACT  has  one 
dominant  caustic  per  convergence  zone.  Since  this  caustic  is  so  nearly  horizontal,  its  illumi¬ 
nated  side  exceeds  ranges  for  which  rays  arrive  from  the  same  family.  Thus,  FACT  exhibits 
sharp  increases  in  loss  at  50  and  97  kyd  at  the  end  of  each  convergence  zone. 

The  widths  of  PLRAY’s  caustics  are  narrower  than  LIRA’S  and  correspond  more 
closely  to  the  widths  of  the  normal  mode  caustic.  However,  the  positions  of  the  peaks  of 
PLRAY’s  caustics  are  compressed  in  the  range  scale  compared  to  the  normal  mode  peaks. 
The  positions  of  the  caustic  peaks  in  LIRA  are  closer  to  the  normal  mode  predictions. 

The  widths  of  the  caustics  increase  from  first  to  second  convergence  zones  for  both 
LIRA  and  PLRAY;  AP2  does  not  show  this  feature. 

CASE  13. 

SSP  7  is  used  with  source  depth  at  100  ft,  target  depth  at  300  ft,  and  frequency  at 
100  Hz.  Propagation  loss  is  presented  for  LIRA,  PLRAY,  and  FACT  in  Fig.  D-15. 

All  three  models  show  propagation  loss  curves  made  up  of  nearly  horizontal  caustics 
with  cutoffs  on  the  illuminated  side  for  ranges  for  which  there  are  no  ray  arrivals.  The  slow 
rises  with  increasing  ranges  are  due  to  the  caustic  fields  in  the  shadow  zone,  and  the  choppy 
nature  of  the  curves  thereafter  reflects  the  cutoff  ranges  from  the  several  ray  families. 

LIRA  has  additional  peaks  at  36  and  72  kyd,  which  are  due  to  rays  converging  but 
not  crossing  one  another  in  range. 

CASE  14. 

This  is  the  same  as  the  previous  case,  with  the  source  depth  changed  to  1 000  ft. 
Incoherent  propagation  losses  for  LIRA,  PLRAY,  and  FACT  are  presented  in  Fig.  D-16. 
Smoothed  coherent  propagation  loss  for  AP2  is  presented  in  Fig.  D-7c. 

The  caustics  calculated  by  PLRAY  and  FACT  occur  at  the  same  ranges.  The  LIRA 
curve  of  propagation  loss  seems  slightly  expanded  in  comparison.  The  corresponding  AP2 
smoothed  coherent  propagation  loss  curve  places  the  caustics  at  the  same  ranges  as  LIRA. 

LIRA  and  FACT  show  the  direct  path  extending  out  to  at  least  1 5  kyd;  PLRAY 
achieves  only  3  kyd.  There  is  a  caustic  at  3  kyd,  which  LIRA  merges  with  the  other  ray 
contributions.  Unfortunately,  the  tail  of  the  illuminated  side  of  the  caustic  is  not  well 
merged  with  the  corresponding  ray  arrivals  for  that  family,  giving  the  dip  at  5  kyd  and  a 
bump  at  6  kyd.  LIRA’S  bump  at  6  kyd  does  not  correspond  to  FACT’S  small  peak  at  3  kyd, 
which  is  the  correct  position  of  the  caustic.  From  the  shape  of  FACT'S  peak,  it  appears  that 
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Figure  D-14.  Incoherent  propagation  loss,  CASE  12  (NADC  SSP  6,  frequency  50  Hz,  source  depth 
1000  ft,  target  depth  300  ft). 


71 


PROPAGATION  LOSS  (OBI 


Figure  D-I5.  Incoherent  propagation  loss,  CASE  13  (NADC  SSP  7,  frequency  100  Hz,  source  depth 
100  ft,  target  depth  300  ft). 


FACT  did  not  actually  find  a  caustic  at  3  kyd,  but  computed  propagation  loss  by  means  of 
ray  theory. 

PLRAY  shows  a  small  peak  at  32  kyd,  which  should  have  been  computed  using  caus¬ 
tic  corrections. 

The  chop  in  the  LIRA  curve  from  40  to  52  kyd  is  from  rays  that  strike  the  surface. 
Propagation  loss  from  each  of  the  ray  types  looks  like  a  plateau  about  8  kyd  in  width  with 
steep  sides.  The  different  ray  types  have  different  ranges  at  which  onset  of  the  plateau  occurs 
as  well  as  the  dropoff.  When  these  plateau-shaped  curves  are  combined  in  the  incoherent 
summation,  the  result  is  the  choppy  curve. 

CONCLUSIONS 

LIRA  shows  generally  good  agreement  with  PLRAY,  FACT,  and  AP2  in  the  direct- 
path  and  convergence  zone  regions  and  excellent  agreement  in  bottom  bounce  regions. 

Where  serious  differences  occur  LIRA  agrees  well  with  at  least  one  of  the  other  three  models. 

LIRA  is  able  to  resolve  more  ray  theory  caustics  than  PLRAY  or  FACT.  FACT 
appears  to  merge  ray  types  to  produce  only  one  caustic  where  two  should  exist.  PLRAY’s 
caustic  peaks  sometimes  appear  steep  and  pointed  and  the  range  scale  compressed.  On  the 
other  hand,  FACT  seems  to  coalesce  two  or  more  caustics  into  a  single  broadened  peak, 
with  the  range  scale  expanded.  The  positions  of  the  peaks  in  LIRA  and  AP2  fall  between 
these  trends. 

Of  the  three  ray  models,  LIRA  and  PLRAY  are  the  two  which  show  closest  agree¬ 
ment  and  the  most  similarity  in  the  propagation  loss  curves. 

LIRA  exhibited  no  pathological  deviations  except  at  21  kyd  in  CASE  1 1 ;  yet 
this  deviation  is  explained. 

On  the  basis  of  these  14  comparison  cases,  LIRA’S  propagation  loss  program  should 
be  accepted  as  a  reasonable  predictor  suitable  to  be  used  for  active-sonar  performance 
prediction  at  low  frequencies. 


APPENDIX  E 


EXAMPLES  OF  OUTPUT  PLOTS  FROM  LIRA 

The  following  plots  demonstrate  the  use  of  LIRA  to  analyze  an  active-surveillance 
scenario.  The  winter  and  summer  sound-speed  profiles  in  Figs.  E-1  and  E-2  are  representative 
of  an  area  at  the  edge  of  a  continental  shelf,  where  the  bottom  drops  sharply  from  3600  to 
17600  ft.  An  active  source  at  700  ft  is  positioned  above  a  horizontal  receiving  array  on  the 
bottom  at  3600  ft. 

Four  cases  are  investigated: 

CASE  SSP  TARGET  DEPTH  (FT.) 

E-1  Winter  60 

E-2  Winter  400 

E-3  Summer  60 

E-4  Summer  400 

Three  plots  are  produced  for  each  case:  propagation  loss,  pulse-averaged  reverberation, 
and  signal  excess. 

The  input  card  images  used  to  run  the  LIRA  program  are  given  in  Table  E-1,  an 
example  of  one  execution  of  LIRA  containing  four  runs.  Input  parameters  of  interest  are: 

•  Frequency  is  250  Hz  (FREQ) 

•  Surface  loss  per  bounce  is  0.0 1 5  dB  (HS) 

•  Volume  column  backscattering  strength  is  -70  dB/yd"//l  yd  (MUV) 

•  Surface  backscattering  strength  is  a  function  of  grazing  angle:  for  low  angles  it 
is  -60  dB/yd“//l  yd  (ASBS). 

•  Bottom  backscattering  strength  is  a  function  of  grazing  angle;  for  low  angles 
it  is -35  dB/yd2//  1  yd  (ABBS). 

•  Source  level  is  240  dB  //  1  pPa  at  1  yd  (SL). 

•  Source  vertical  beam  pattern  is  that  of  a  continuous  line  array  (BEAMX  =  0). 

•  Vertical  beamwidth  of  the  main  lobe  of  the  source  beam  pattern  is  20  deg  (VBM). 

•  The  signal  waveform  is  PRN. 

•  The  pulse  length  is  10  s  (PULSE). 

•  Receiving  array  gain  is  1 3  dB  (AG). 

•  Receiving  array  vertical  beam  pattern  is  omnidirectional. 

•  Receiver  bandwidth  is  100  Hz  (BWR). 

•  Spectrum  level  noise  is  64  dB/ 1  Hz  //  1  pPa  (NIN). 

•  Target  strength  is  1  5  dB  (TGS). 

•  Detection  threshold  is  14.5  dB  (ADT)  for  0.9  probability  of  detection  and 
10"4  probability  of  false  alarm. 

•  Only  up-rays  are  allowed  at  the  receiver  (RBOT). 

Propagation  loss  vs  target  range  for  the  four  cases  is  presented  in  Figs.  E-3  to  E-6. 

The  propagation  loss  shown  is  the  average  of  the  one-way  losses  from  source  to  target  and 
from  target  to  receiving  array.  The  solid  curve  is  the  incoherent  sum  of  all  ray  arrivals,  and 
the  dashed  curve  is  the  propagation  loss  for  the  largest  arrival  (single  path  with  least  loss). 

The  beam-pattern  effects  are  included  in  the  loss  curves. 
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LIRA 


LOW-FREQUENCY  INTERMEDIATE-RANGE  ACTIVE-SONAR  PERFORMANCE  PREDICTIONS 


RESPONSIBLE  PERSON:  DAVID  *.  HOFFMAN 

NAVAL  OCEAN  SYSTEMS  CENTER 
CODE  7142 

SAN  DIEGO,  CA  92152 

TELEPHONE  NUMBERS  (714)  225-2316 

AUTOVON  935-2316 


•  INPUTS,  RUN  1  * 


HEADER 

CASE  E -  1 ,  WINTER 

SSP  0  1499  bO  1501  140  1493  200  1490.5  350  14b5  4o0  U82.5  600  1480 

•  995  1482  1  500  1487.5  2000  U93  3000  1  5 Oo . 5  5200  1549 

FREQ  25C.SL  240, PULS  10.TGS  15 

I C  1000, MUV  -70 
N1N  67, DI  13.BWR  100 
HS  *  1  .WAVE  3. 5, oINO  10 
ADT  14.5 

AHB  0  3  6  3  e  3.14  10  3.71  12  4.71  14  6.13  16  7.56  18  8.98 

•  20  10.4  22  11.32  24  12.24  26  12.9  26  13.3  3u  13.7  32  13.85  34  14 

•  90  14 

AS8S  0  -60  10  -60  30  -40  60  -30  80  0  90  5 

ABBS  0  -35  77  -27  80  -15  65  -2  90  0 

8  E  AM  X  O.DELPX  20, X  DE  0 

ORONLY,  LA,  RBOT,  TTY 

PLOT  1  1  1 

R  1  2  600  2 

2  R  3600 

ZX  700, ZTG  60, RUN 


*  CALCULATIONS,  RUN  1  * 


DETECTION  ANNULUS  FROM  .0  TO  5.5  AYD 
DETECTION  ANNULUS  FROM  61.7  TO  63.6  AYD 


Tabic  E-l.  An  example  of  one  execution  of  LIRA  with  four  runs. 


*  INPUTS,  RUN  2  * 


HEADER 

CASE  E-2,  PINTER 
2  TG  AGO , RUN 


*  CALCULATIONS,  RUN  2  * 


DETECTION 

ANNULUS 

FROM 

,  G 

TO 

5.8 

K  YD 

DETECTION 

ANNULUS 

FROM 

<.1.8 

TO 

A3. 3 

KYD 

DETECTION 

ANNULUS 

FROM 

<.7.7 

TO 

62.0 

K  YD 

DETECTION 

ANNULUS 

FROM 

99.7 

TO 

101.7 

KYD 

DETECTION 

ANNULUS 

FROM 

103.3 

TO 

107.2 

K  YO 

DETECTION 

ANNULUS 

FROM 

109.9 

TO 

113.5 

KYD 

DETECTION 

ANNULUS 

FROM 

117.9 

TO 

119.2 

KYD 

DETECTION 

ANNULUS 

FROM 

155.1 

TO 

157.9 

K  YO 

DETECTION 

ANNULUS 

FROt 

161.5 

TO 

165.6 

KYD 

DETECTION 

ANNULUS 

FROM 

215.7 

TO 

217. A 

KYD 

DETECTION 

ANNULUS 

FROM 

225.6 

TO 

229.2 

KYD 

Table  E-l.  Continued. 
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*  INPUTS 


RUN  3  * 


HEADER 

CASE  E-3  f  SUMMER 

SSP  0  151a  30  1516.5  50  1505  70  1500  100  1495  170  1488  400  1480.5  500  1479 

*  600  1477.5  1000  1479  1500  U84  2000  1490  3000  1506  5200  1546 
ZX  700, ZTG  60, RUN 

••»*»*»*«*»**«*»** 

*  CALCULATIONS,  RUN  3  « 

***********4************* 


DETECTION 

ANNULUS 

FROM 

.0 

TO 

5.3 

K  YD 

DETECTION 

ANNULUS 

FROM 

55.6 

TO 

57.9 

K  YD 

DETECTION 

ANNULUS 

FROM 

59.1 

TO 

61  .3 

K  YD 

DETECTION 

ANNULUS 

FROM 

63.7 

TO 

65.7 

KYD 

DETECTION 

ANNULUS 

FROM 

113.6 

TO 

115.4 

KYD 

DETECTION 

ANNULUS 

FROM 

117.7 

TO 

119.2 

KYD 

DETECTION 

ANNULUS 

FROM 

121.7 

TO 

123.4 

KYD 

Table  E-l.  Continued. 


*  INPUTS,  RUN  4  * 


HEADER 

CASE  E-4,  SUMMER 
Z  T G  400, RUN 


*  CALCULATIONS,  RUN  4  * 
******  ***** ******  ******** 


DETECTION 

ANNULUS 

FROM 

.0 

TO 

5.2 

KYD 

DETECTION 

ANNULUS 

F  ROM 

41 .9 

TO 

43.3 

K  YO 

DETECTION 

ANNULUS 

FROM 

47,6 

TO 

59.6 

KYD 

DETECTION 

ANNULUS 

FROM 

99.5 

TO 

107.9 

KYD 

DETECTION 

ANNULUS 

FROM 

1G9.0 

TO 

113.7 

KYD 

DETECTION 

ANNULUS 

FROM 

115.5 

TO 

117.3 

KYD 

DETECTION 

ANNULUS 

FROM 

149.9 

TO 

158.0 

KYD 

DETECTION 

ANNULUS 

FROM 

159.5 

TO 

164.0 

K  YD 

DETECTION 

ANNULUS 

FROM 

165.9 

TO 

167.1 

KYD 

DETECTION 

ANNULUS 

FROM 

201.6 

TO 

209.2 

KYD 

DETECTION 

ANNULUS 

FROM 

213.1 

TO 

215.5 

KYD 

DETECTION 

ANNULUS 

FROM 

253.4 

TO 

255.8 

KYD 

DETECTION 

ANNULUS 

FROM 

263.8 

TO 

265.7 

KYD 

DETECTION 

ANNULUS 

FROM 

303.9 

TO 

305.1 

KYD 

DETECTION 

ANNULUS 

F  ROM 

307.9 

TO 

30V. 1 

KYD 

DETECTION 

ANNULUS 

FROM 

410.0 

TO 

411.1 

KYD 

*  INPUTS, 

*********1 

k  *******  1 

RUN  5 

k* 

* 

k  * 

END 


Table  E-l.  Continued. 


81 


PROPAGATION  LOSS  VS  RANGE 

OUTPUT  PROM  LIRA 
CASE  E-l,  WINTER 


Figure  E-3.  Propagation  loss  for  CASE  E-l 
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OUTPUT  FROM  LIRA 
CASE  E-4,  SUMMER 
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Reverberation  level  vs  target  range  for  the  four  cases  is  presented  in  Figs.  E-7  to 
E-10.  The  reverberation  level  corresponds  to  the  time  the  leading  edge  of  the  signal  arrives 
at  the  receiving  array  after  bouncing  off  the  target  at  the  specified  range.  The  plotted  rever¬ 
beration  level  in  dB  //  1  ^Pa  is  averaged  over  pulse  length  and  corrected  for  processing  gain. 

This  corresponds  to  reverberation  in  the  receiver  bandwidth  at  the  output  of  the  signal 
processor  before  the  thresholding  device.  Each  of  the  bottom,  surface,  and  volume  com¬ 
ponents  of  reverberation  is  the  incoherent  sum  of  the  respective  intensities  derived  from 
ray  tracing. 

Noise  threshold  displayed  on  the  reverberation  plots  is  the  ambient  noise  and  self¬ 
noise  in  the  receiver  bandwidth  at  the  output  of  the  signal  processor.  In  this  example  noise 
threshold  is  calculated  by  LIRA  by  correcting  spectrum  level  noise  for  array  gain,  band¬ 
width,  and  processing  gain. 

Signal  excess  vs  target  range  for  the  four  cases  is  presented  in  Figs.  E-l  1  to  E-14.  A 
signal  excess  of  zero  corresponds  to  0.9  probability  of  detection  with  1 O'4  probability  of 
false  alarm.  The  solid  curve  is  the  signal  excess  for  propagation  loss  to  the  target  made  up 
of  the  incoherent  sum  of  all  ray  arrivals.  This  curve  represents  optimistic  performance  ob¬ 
tainable  only  if  all  the  target-return  echoes  for  the  multipaths  can  be  resolved  and  recombined. 
The  solid  curve  represents  an  upper  limit  on  sonar  performance.  The  dashed  curve  is  the  proper 
signal  excess  curve  to  use  for  sonar  performance.  It  corresponds  to  the  largest  arrival,  the 
strongest  echo  on  any  single  path. 

LIRA  has  the  option  to  print  detection  ranges,  i.e.,  annulae  in  which  the  signal  excess 
exceeds  zero.  Detection  ranges  for  the  four  cases  are  given  in  Table  E-l . 


AVERAGED  REVERBERATION  VS  RANGE 

OUTPUT  FROM  LIRA 
CASE  E-2,  WINTER 


SIGNAL  EXCESS  VS  RANGE 

OUTPUT  FROM  LIRA 
CASE  E-l,  WINTER 
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